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ABSTRACT 

Aims. We present 139 spectra of 124 Type la supernovae (SNe la) that were observed at the ESO/VLT during the first three years 
of the Canada-France-Hawal Telescope (CFHT) Supernova Legacy Survey (SNLS). This homogeneous data set is used to test for 
redshift evolution of SN la spectra, and will be used in the SNLS 3rd year cosmological analyses. 

Methods. Spectra have been reduced and extracted with a dedicated pipeline that uses photometric information from deep CFHT 
Legacy Survey (CFHT-LS) reference images to trace, at sub-pixel accuracy, the position of the supernova on the spectrogram as 
a function of wavelength. It also separates the supernova and its host light in ~60% of cases. The identification of the supernova 
candidates is performed using a spectrophotometric SN la model. 

Results. A total of 124 SNe la, roughly 50% of the overall SNLS spectroscopic sample, have been identified using the ESO/VLT 
during the first three years of the survey. Their redshifts range from z = 0.149 to z = 1.031. The average redshift of the sample is 
z = 0.63 ± 0.02. This constitutes the largest SN la spectral set to date in this redshift range. The spectra are presented along with their 
best-fit spectral SN la model and a host model where relevant. In the latter case, a host subtracted spectrum is also presented. We 
produce average spectra for pre-maximum, maximum and post-maximum epochs for both z < 0.5 and z > 0.5 SNe la. We find that 
z < 0.5 spectra have deeper intermediate mass element absorptions than z > 0.5 spectra. The differences with redshift are consistent 
with the selection of brighter and bluer supernovae at higher redshift. 

Key words. cosmology:observations ~ supernovaeigeneral ~ methods: data analysis ~ techniques: spectroscopic 



1. Introduction 

Since the direct detec tion of the acceler ated expansion of th e 
universe 10 years ago (iRiess et al . 1998; P erlmutter et al.|[r999l) . 
constraining the equation of state of the dark energy compo- 
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* Based on observations obtained with FORSl and FORS2 at the 
Very Large Telescope on Cerro Paranal, operated by the European 
Southern Observatory, Chile (ESO Large Programs 171.A-0486 and 
176.A-0589) 

** Figures A.l to A. 139 are only available in electronic form via 
|http://www.edpsciences.org| 



nent responsible for this acceleration has been a major goal 
of observational cosmology. Type la supernovae (SNe la here- 
after) samples have been gathered at low and high redshift and 
extensively used for this purpose. When combined with other 



emerges (iTonrv et al.ll2003l; lAstier et al.l)2006 


IRiess et al.ll2007l: 


IWood-Vasev et al.ll2007UKowalski et al.ll2008 


). 



Over the past five yearfl the Supernova Legacy Survey 
(SNLS) has gathered more than 1000 light curves of SN la can- 
didates on the Canada-France-Hawai telescope (CFHT) using 



The SNLS started in June 2003. 
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MEGACAM dBoulade et al.ll2003l) . thanks to a rolling search tech- 
nique for discovery and photometric follow up of SNe la in four 
1 square degree fields (Astier et al. 2006). Spectra of a subset 
(about half) of these SNe la candidates have been observed on 
various 8 to 10-m class telescopes ( VLT, Gemini N and S, Keck I 
and II). About 50% of spectroscopically observed SN candidates 
were observed at the VLT. 

In this paper, we present the VLT SN la spectral set for the 
first three years of operation of the SNLS. The non SN la spec- 
tral set, together with a description of the "rea l-time" operation s 
and procedures will be presented elsewhere (iBasa et al.ll2009l) . 
The spectra shown here were obtained in the period running 
from June 1st, 2003 up to July 31st, 2006, as part of two large 
VLT programs Q. The spectra have been analysed using novel 
techniques for extraction and identification which have been 
described in detail elsewhere (i Baumonj 120071 : iBaumont et al.l 
|2008|) . For each spectrum we provide a redshift estimate. The 
identification of SNe la relies in pa rt on huma n judgement, us- 
ing the SALT2 spectral te mplate of iGuv et alj ( 12007.) as a guide 
(see IBaumont et al.ll2008h . 

The SN spectra presented here will be used, along with SN la 
spectra obtained at Gemini and Keck telescopes, to build the 3rd 
year SNLS Hubble diagram. Our primary goal is that the type 
and redshift of the SNe la used for cosmological analysis are 
secure. In this paper, we consider two classes of events: secure 
SNe la ("SN la") and probable SNe la ("SN la*"). Studying the 
statistical properties of these two classes, in order to assess the 
validity of using SN la* events together with SN la events, is 
therefore a goal of this paper. 

SN la spectra are a rich source of physical information 
about their progenitor history and environment. The possibil- 
ity of evolution among SN la populations at low and high red- 
shifts has been the subject of considerable attention in recent 
years, as more and m ore SN la spec tral sets become available 
fcaravini et al.' '2007": 'Bronder et al.' '20081 iFolev et alj l2008at 
[Blondin et al. 2006; Ellis et al. ^008). Recently, evidence has 
even been found for a demographic evolution among SN la pop- 
ulations, resulti ng in higher stretch, more luminous SNe la at 
higher redshift jHowell et a l. 2007; Sullivan et al. 2009). Using 
SN la spectra to compare their physical properties at low- and 
high-redshift is therefore a useful cross-check when using SN la 
to constrain the expansion history of the universe. In this pa- 
per, we take advantage of the large number of high quality spec- 
tra obtained at the VLT to build average composite spectra at 
various phases with respect to maximum light for z < 0.5 and 
z > 0.5. We also compare our average spectra to comp osite spec- 
tra obtained in a similar way with different data sets dEllis et alj 
120081: [Foley et al.l l2008a') and discuss the significance of the dif- 
ferences found in terms of possible evolution or selection effects. 

A plan of the paper follows. In Section|2l we briefly describe 
the SNLS photometric survey and the VLT spectral observation 
programs. In Section[3j we summarise the main steps of the data 
reduction and spectrum extraction. We detail our identification 
procedure and classification scheme in Section ID In Section |5] 
the spectra are individually presented. Composite spectra at z < 
0.5 and z > 0.5 are built in Section |6] In Section |7] we discuss 
our sample in the light of other existing SN la spectroscopic data. 
Concluding remarks are made in Section[8] 



^ These dates correspond to the first three years of the SNLS. 



2. SNLS Observations 

2.1. The SNLS imaging survey 

The SNLS is composed of an imaging survey devoted to the de- 
tection and the photometric follow up of SN candidates, and 
a spectroscopic survey of a sample of the detected candidates, 
prioritised for spectroscopy on various telescopes. The imaging 
survey ran from June 2003, after a period of pre-survey, until 
June 2008. It was based on the Deep survey of the Canada- 
France-Hawaii Telescope Legacy Survey (CFHT-LS) (amount- 
ing to half of the 474 nights allocated to t he CFHT-LS). FuU 
details of the survey can be found elsewhere d Astier et al.ll2006l) . 
In brief, SNLS observed 4 fields (Dl-4) every 3-5 nights during 
dark/grey time in the griz filters, each field followed for 5-6 lu- 
nations per year. Around 1000 well sampled multi-colour light 
curves of SN la candidates have been obtained up to z ~ 1.1. 

2.2. Spectroscopic follow-up 

Spectroscopy of SNLS SN candidates was performed on several 
8 to 10-m class telescopes in both hemispheres, namely the VLT, 
Gemini-N and S, Keck I and II. Almost 50% of SNLS candidates 
identified as certain or probable SN la we re spe ctroscopically 
observ ed at the VLT. iHoweU et all ( l2005h and iBronder et al] 
( |2008|) describe the SNLS fi rst three years of Gemini spectral 
data (up to May 2006), while lElfis et alj (l2008l) present 36 high 
signal-to-noise ratio (S/N) SN la spectra obtained at Keck. In this 
paper, we focus on spectra taken at the VLT on Cerro Paranal. 

Candidate selection for spectroscopic follow -up was based 
on the multi-band photometry procedure of ISullivan et al.l 
(2006). This selection was performed on the rising part of the 
light curves, routinely available thanks to the rolling search strat- 
egy (see e.g.. |Perrett et al.l2009b . Candidates were generally sent 
for spectroscopy at, or slightly after, maximum light, which op- 
timised the time budget allocated for spectroscopy. We triggered 
a target on VLT every three to four days during dark and grey 
time. 

During the first large program (2003-2005), we performed 
long slit spectroscopy (LSS) of SN candidates on FORSl for a 
total of 60 hours of dark/grey time per semester. During the sec- 
ond large program (2005-2007), we observed using both FORS 1 
and FORS2 with the standard collimator in LSS and multi object 
spectroscopy (MOS) mod^ Most observations were carried out 
with the 300V grism, along with the GG435 order-sorting filter 
Grism 300V was chosen to optimise spectral resolution, spectral 
coverage and high enough S/N for an unambiguous identifica- 
tion, even for the faintest candidates of our survey (i'^g ~ 24). 
Moreover, using the 300V grism for high redshift SNe allows us 
to study the interesting rest frame UV region of the spectra. The 
pixel scale is 0.2" along the spatial axis and 2.65 A along the 
dispersion axis. At 5000 A, the resolution limit with this setup 
is 1 1 A. The efficiency of the 300V grism peaks around 4700 
A at a level of ^ 80%. The 3001 grism along with the OG590 
order-sorting filter was sometimes used for the faintest (z > 0.8) 
SNe. Being typically more distant, and due to strong sky emis- 
sion and fringing beyond 6000 A, spectra obtained in this way 
have a much lower S/N than those acquired with the 300V grism. 



^ Only eight candidates from SNLS3 were observed in MOS mode 
and have been identified as SN la in real-time. As the MOS mod e is cur- 
rently not supported by our new extraction pipeline (see Baumont et alj 
120081) . we do not include them in here. Note also that only 3 SNe la 
were observed on F0RS2 in the period covered in this study. 
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The slit width was chosen according to the following rule 
of thumb: "slit width » seeing + 0.2"", as a compromise be- 
tween observing most of the flux from the targeted candidate and 
limiting the sky background flux. An air mass < 1 .6 was re- 
quired for each spectroscopic observation. A "blind offset" tech- 
nique was used to target the candidate, using a bright star located 
within < 1 ' of the target and then offsetting the telescope to po- 
sition the slit onto the candidate. When possible, the slit position 
was chosen to observe both the SN and its host. Differential slit 
losses were corrected by a Longitudinal Atmospheric Dispersion 
Corrector (LADC). Residual losses are taken into account with 
the recalibration procedure described in Section l4!2l 

All spectra were acquired in Service Observing mode. With 
a limiting magnitude of i'^^ ~ 24, 3-4 exposures of 750 or 900s 
were taken for each candidate, with small offsets along the spa- 
tial axis (Y; the dispersion axis X is horizontal with our setup). 
Thanks to the regular time sampling of the rolling search, it was 
possible to acquire most of candidates around or slightly past 
maximum light. 

3. Data processing 

Data reduction and spectral extraction were performed in two 
separate ways. A quick, "real-time" analysis (within a day of 
acquisiti on) was used to assess the type and redshift of the can- 
didate ( Basa et aLll2009l) . an essential task for efficiently allo- 
cating other candidates to the various telescopes. In parallel, 
we developed tools "ofi'-line" to cleanly extract the SN from 
its host. A dedicated pipeline, PHASE (PHotometry Assisted 
Spec tral Extraction), was used for the final reduction and extrac- 
tion (lBaumontll2007l: iBaumont et all 12008) . Full details can be 
found in these papers; we give only a brief description here. All 
extractions presented in this paper used the PHASE technique. 

The PHASE reduction technique improves over the real-time 
reductions in refinements of the master flat-fields, the dispersion 
relation, and in the sky estimation. As an example, the 2D disper- 
sion relation is modeled by a fourth order polynomial in X, (and 
of 2nd order in F, and XYi) to further reduce the residuals. For 
flux calibration, we build a single average response curve (one 
for UTl and one for UT2 as FORS 1 moved from UTl to UT2 in 
June 2005) from previous individual standard star observations. 
We prefer using a well controlled average response for the whole 
set, rather than using a response built from standard star obser- 
vations of a different night. This is at the expense of absolute 
flux calibration, as we average out sky transmission variations 
from night to night, but permits a more robust estimation of the 
sensitivity function near the blue edge of the order sorting filter. 

PHASE extraction uses photometric information on the SN 
location and host brightness from the deep reference images 
(/^g ~ 27 - 28) used for building the SN light curvefl This 
allows us to accurately trace the SN position along the disper- 
sion axis on the two-dimensional spectrogram. Moreover, we 
build a multi-component model of the galaxies present in the slit 
(including the host, if resolved) by measuring the spatial pho- 
tometric profiles of these galaxies on the deep stacked reference 
images projected along the slit direction. We then add a SN com- 
ponent, modeled as a Gaussian of width equal to the seeing of the 
spectroscopic observation. The location of the SN is accurately 
known from the light curves. The flux of each component (SN + 
host and potentially other galaxies in the slit) is a free parameter, 
the sum of the profile fluxes being normalised to unity. Such a 

We use gM, rM and z'm for spectra obtained with grism 300V 
dBaumont et a].ll2008h 



model is built for each column (the dispersion direction is hori- 
zontal with our setup) and fluxes assigned to each component in 
each column / are estimated by ax^ minimisation where; 

j 

Here, o-j is the error associated with the column pixel data, Aii 
is the multi-component model in column / and Fi{Yj) = F(X,, Y) 
is the column flux on the 2D spectrogram. We distinguish three 
types of galaxies: 

- PSF: unresolved, point-Uke galaxies 

- EXT: extended, but regularly shaped profiles, e.g. ellipticals 

- Mix: extended, but irregularly shaped profiles, e.g. galaxies 
with spiral arms. 

We note that the latter two cases are not an accurate indica- 
tion of the physical morphology of the host galaxy, but instead 
just model the spatial profile of the host that enters the slit. 

In the first case, the galactic component of the model is a 
Gaussian of width equal to the seeing of the deep reference 
observation. The seeing variation with wavelength is estimated 
from standard star spe ctra as a power law of index ~ -0.3, 
in good agreement with iBlondin et al.1 (l2005b measurements on 
FORS 1 spectra. In the second case, the model is the "bolomet- 
ric" spatial profile (the sum of the galactic profile in all observed 
filters) as measured on the megacam combined deep reference 
images. For the third case, we use a mixture of a Gaussian with 
width equal to the spectroscopic seeing to model the core and the 
photometric spatial profile to model the extended arms. From a 
pure algorithmic point of view, this latter case is equivalent to 
have two distinct galaxies, a point-like source and an extended 
one. 

Host models used to extract our VLT spectra are about 
equally divided into EXT and PSF models (~ 30% each), with 
only a few percent of cases being Mix models. In the remaining 
~ 40% of cases, no separate extraction of the SN and the host 
was possible. 

PHASE uses a set of default parameters to select the correct 
host model and make the extraction as automatic as possible. 
These parameters include cuts on flux, galactic compactness, 
extension minimum level, colour variation between the centre 
and the extended part of the host (identifying "Mix" host types), 
and a minimum value of the SN to host centre distance to per- 
form a separate extraction (usually 0.15", a bit less than one 
pixel). These default parameters allow an automatic extraction 
of most spectra, though the parameters can be adjusted for spe- 
cific case s. PHASE performanc e and limitations have been dis- 
cussed in iBaumont et al.l (l2008h and we refer the reader to that 
paper. 

The main hypotheses in using PHASE are that 1) the PSF 
is a Gaussian of width equal to the seeing; 2) the coordinates 
of the SN are accurate; 3) CFHT-LS reference images and VLT 
spectrograms have comparable seeings; 4) no flux of the SN is 
present in the reference image. Any deviation from these as- 
sumptions result in weak flux losses, increased noise, and con- 
tamination of a SN spectrum by its host. PHASE performs well 
for most spectra encountered. In particular, it succeeds at recov- 
ering the SN from the host, even in the case of a SN located 
close to the host centre (typically J 1 pixel): strong correlations 
between the host and the SN are often unavoidable in standard 
extractions. Even for "non favourable" cases, such as sub-pixel 
SN/host separations, both component spectra are recovered and 
are essentially non correlated. This is a remarkable feature of our 
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pipeline, as most SN spectral extractions are hampered by host 
contamination in these cases. Nevertheless, if the SN is too close 
to the host centre (separation less than one fifth of the seeing), 
no separate extraction is possible. In that case, a host spectral 
template is used at the stage of identification to estimate the host 
contamination (see Section l43T l. 

A comparison between PHASE and standard extracti ons has 
been d one and illustrated on a few examples in Baum ont et alj 
(l2008h . A major difference is that, in PHASE, we do not re- 
sample the data until the final step of flux calibration. This avoids 
introducing correlations across pixels and allows us to trace the 
statistical noise along the reduction and extraction procedure. 
For the same reason, we also avoid rebinning the data at a con- 
stant wavelength step, as is done in most standard spectroscopic 
reduction pipelines. As a consequence, the final PHASE cali- 
brated spectrum has unequal steps. We have checked that the 
statistical noise is properly propagated with PHASE by comput- 
ing < S /N >, the S/N per pixel averaged over the whole spectral 
range, and < S /cr >, the signal-to-r.m.s. ratio. Here, cr is the 
standard deviation of a group of measurements around a low- 
order polynomial fit. We find relatively good agreement between 
these two quantities for PHASE extracted spectra. 

4. Spectral analysis 

4.1. Redshift determination 

Where possible, a redshift is obtained from the host galaxy spec- 
trum using characteristic emission or absorption lines, yielding 
a typ ical uncertainty of 0.001 on the redshift (iBaumont et al.l 
12008) . When no lines are detected, the redshift is determined 
from the fit of a model to the SN spectrum (see Section |4.2| i. 
The typical redshift uncertainty is then Sz ~ 0. 01, due to the di- 
versi ty of ejecta velocities among SNe la, (e.g. iHachinger et al.l 
12006 ). In alDout 80% of cases, the redshift is obtained from host 
emission and/or absorption lines. In SNLS we have two indepen- 
dent identification pipelines. The redshifts (as well as identifica- 
tions, see below) have been carefully cross-checked using the 
two-dimensional data to match host lines. Checking the corre- 
sponding noise map, "bright" spots visible on the 2D frames are 
easily identified as true emission lines or cosmic/sky subtraction 
residuals. 

4.2. Identifying SNe la 

To identify SNe la, we use a minimisa t ion pr ocedure using the 
SALT2 spectral template of iGuv et al.l (|2007[) with a combined 
fit of the hght curves and the spectrum consistently performed 
dBaumont et all l2008l) . As the training sample of the SALT2 
model only contains SN la spectra and light-curves, it does not 
allow a direct identification of non-SN la objects. However, the 
best-fitting parameter values obtained when fitting a non-SN la 
with SALT2 are in themselves an indirect indication of the SN 
tvpe (see Section |5.4| i. In addition, a template-fitting code 
(iHowell et alJl2005l) has been used for cross-check. Both tech- 
niques use varying levels of human judgement. 

The main output parameters of the SALT2 fit are 1) the light 
curve fit parameters xq (overall normalisation), xi (light curve 
shape), and colour c, and 2) spectroscopic fit parameters: the host 
fraction fgai in the model when relevant, and recalibration pa- 
rameters. The latter enter a recalibration function applied to the 
photometric model in order to fit the spectru m and account for 
possible errors in flux calibration (iGuv et alj|2007l) . This func- 
tion is a polynomial of order n, with coefficients y„ inside an 



exponential (to ensure positivity). We usually only use two re- 
calibration parameters: an overall normalisation yo and a first 
order coefficient yi (tilt applied at rest frame 4400 A), adding 
higher order corrections for a few objects. The xi parameter is a 
light curve shape parameter that can be c onverted int o a stretch 
factor s or Amis parameter of iPhillipsI (|T993); see iGuv et all 
(2007). The colour c parameter is defined as the difference be- 
tween (B - V)s N and the average < B -V > value at maximum 
light for the whole training sample of SALT2. 

There are some advantages in using SALT2 as a tool for 
identification. First, the fit of the light curve yields the date of 
maximum light. The phase (p of the spectrum (the number of 
rest frame days between the date of B-band maximum and the 
date of acquisition of the spectrum) is accurately known, usu- 
ally within a fraction of day. This alleviates possible degenera- 
cies between SN types, such as between pre-maximum Type Ic 
supernovae (SNe Ic) and post-maximum SNe la, whose spectra 
show similarities. Second, the training set of the SALT2 model 
is built from a large collection of spectra and light curves from 
local and distant SNe la. These latter include SNLS SNe them- 
selves, added to the training sample once identified. Third, using 
a model instead of a set of spectrum templates (as is typical of 
standard identification techniques) alleviates the problems due to 
unavoidable incompleteness in phase and wavelength coverage 
of the template libraries: all SNe are treated on an equal footing. 
Note also that only SNe typed as SN la (as opposed to SNIa*) 
are included in the training set of SALT2. As an erroneous typ- 
ing of a SN Ic as a SN la is very unlikely, this limits the chance of 
polluting the training sample. Even in such unlikely case, SALT2 
being a model, it is robust against the inclusion of a non-SN la 
object. 

Once the SALT2 fit is performed, the identification is guided 
by the best-fitting spectral parameter values. As the phase is 
fixed by the light curve fit, the spectral fit implicitly uses the 
photometry. However, we will not classify a SN candidate as a 
certain SN la if we do not get an adequate fit of the spectrum, 
even if the light curve fit is good. On the contrary, a convincing 
spectral fit is sufficient for an SN la identification, even if the 
light curve fit is poor. Our goal is to obtain a clean, spectroscop- 
ically confirmed, SN la sample. 

4.3. Host galaxy subtraction 

Spectroscopic identification of distant SNe is challenging. One 
of the key issues discussed in the literature is host contamina- 
tion. Low S/N spectra are common, either due to high redshift 
or because sometimes candidates are observed at a late phase 
(up to a few weeks past maximum) due to telescope schedul- 
ing. Several techniques have been developed to improve host- 
SN s eparation. Standard techniques involve template fitting 
(e.g.. lHowell et al.ll2005l : lLidman et alj|2005l) once the spectrum 
has been extracted, and/or cross-correlation methods, suc h as the 
SuperNova IDentification (SNID) algorithm (Blondin & Tonryl 
20( )7|), used for the first two years of t h e ESSENCE projec t 
CMatheson et al.ll2005HFolev et al.ll2008bh . lBlondin etsU (l2005l) 
proposed a PSF deconvolution technique that separates the two 
components, provided that the spatial extension of the Gaussian 
profile is very different for the SN and for its host. IZheng et al.l 
(llOOSh use a principal component analysis decomposition, com- 
bined with template fitting to assess the level of host contamina- 
tion. 

SNLS has a key advantage in that deep photometric ugriz 
data are available for the host galaxy that can assist with SN-host 
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separation. iBronder et alj (l2008l) used deep iM -band photometry 
to estimate host contamination. More recently. ' Ellis et aLl(l2008l) 
used the multi-colour SNLS photometry to do a more sophisti- 
cated estimate of host contamination in extracting Keck SNe la 
spectra. This approach is very efficient at removing the host con- 
tribution. Though a wide range of host templates are used and are 
colour-matched to the host galaxy photometric data, it neverthe- 
less has the drawback of using synthetic templates. By contrast, 
PHASE, while relying on photometr ic priors, does not use fixed 
host templates dBaumont et al.ll2008h . 

PHASE measures the photometric profile of the host and, 
once combined with a PSF model for the SN placed at the SN po- 
sition, estimates the flux of each component in each pixel along 
the SN trace on the spectrogram using a minimisation proce- 
dure. Inspection of the residual image after extraction shows that 
this technique is efficient provided the centre of the two compo- 
nents is separated by more than ~ 0.15" and th e seeings of the 
reference and spectroscopic images are similar. iBaumont et al.l 
(120081) have shown that fitting a SN spectrum extracted sepa- 
rately from its host, with a SN+host SALT2 model, yields a 
small amount of residual galactic contribution, i.e. less than 10% 
on average. This is a good hint that our component separation 
performs well. Note that our technique is designed within the 
framework of a rolling search (very deep reference images are 
required) and may not be well suited for other search techniques. 

When it is not possible to extract the two components sep- 
arately, ax^ fit of the full (SN+host) spectrum using SALT2 is 
done. SALT2 has been adapted to fit a galaxy template in ad- 
dition to the SN model whenever a separate extraction of the 
SN fr om its host was not possible with PHASE (Baumont et aD 
12008 ). These galaxy templates include iKinnev et ah 1 19961) 
types, as well a s a series of template spectra synthesis ed us- 
ing PEGASE2 ( Fioc & Rocca-Volmerang"3ll997l [1999^ rang- 
ing from ellipticals to late-type spirals. Templates are ordered 
in a regular sequence from red to blue spectra and the best- 
fit model is interpolated between two contiguous templates in 
the sequence. We do not add emission lines to the PEGASE 
templates: whenever PEGASE templates are used to model the 
host contribution, residual emission lines might be found in the 
host subtracted spectrum. This technique is then ess entially com- 
parabl e to the PCA+x^ fitting technique used by I Zheng et al] 
(I2008h to evaluate the host contamination in the first season of 
the SDSS-II survey. 

At high redshift and/or for SNe close to their host centre, 
host subtraction remains a difficulty. We find a higher average 
host fraction for probable SNe la (SN la* type) than for certain 
SNe la (SN la type) spectra (see Section |5T4] l, which shows that, 
as expected, significant host contamination can alter the quality 
of the identification. A major improvement would be to have a 
spectrum of the host. We are currently in the process of obtaining 
"SN free" spectra for those cases in which the subtraction of the 
host failed. We plan to use those to improve the efficiency of the 
host subtraction in the final 5-year SNLS sample. 

4.4. SN Classification 

We classify candid ate spectra into six c ategories adapted from 
the classification of iLidman et al.l (l2005h . We define SN la (cer- 
tain SNe la), SN la* (probable SNe la but other types, in partic- 
ular SNe Ic, can not be excluded given the S/N or the phase of 
the spectrum), SN Ia_pec (peculiar SNe la), SN? (possible SN of 
unclear type), SN Ib/c and SN II. In this paper, we only present 
SNe from the first 3 categories: 



• A SN la classification requires the presence of at least one 
of the following features: Si n /I 4000 or Si ii A 6150, or the S ii 
W-shaped feature around 5600 A. Where these features are not 
clearly visible, a SN la classification is still possible provided 
that the following criteria are met: 

- 1) The overall fit is visually good over the entire spectral 
range, 

- 2) The spectral phase is earlier than about +5 days. At about 
one week past maximum, it has been noted that SNe la and 
SNe Ic show str ong similarities and confusion between types 
is possible (e.g.. lHook et al.ll2005l: iHowell et al.ll2005h . 

- 3) No strong recalibration is necessary to obtain a good fit 
(typical flux correction less than 20%). A large recalibration 
usually indicates that the candidate is not a SN la, a fact that 
would also be reflected in unusual photometric parameters 
(very red or blue colour, i.e. positive or negative c respec- 
tively, very high or low xi value). 

• Classifying a candidate spectrum as a SN la* implies that 
no typical SN la absorption (Si or S) can be found but that the 
overall fit is acceptable over a large spectral range and broad 
features are well reproduced. Spectra of low S/N, or spectra one 
week (or more) past maximum fall into the SN la* category 
unless Si ii is clearly seen. 

• A spectrum is classified as a SN Ia_pec when spectral 
features characteristic of under- or over-luminous objects (e.g. 
iLi et al MoOl) are present. 

While the SN la category represents two thirds of our to- 
tal sample, the number of SNe identified as "peculiar" is very 
small, which may re flect selection biases against such events 
dBronder et al.l l2008h . Low-stretch, underluminous SNe la are 
very hard to identify spectroscopically because of their dim- 
ness compared to their usually bright host. Moreover, low S/N 
spectra of peculiar SNe la at high redshift may make them 
harder to identify than their low-redshift counterparts. Finally, 
our SALT2 -based identifi cation is by essence not well suited for 
handling peculiar SNe la jGuv et al. 2007). 

As described above, all types presented in this paper 
have been cross-checked independently using tw o techniques: 
PHASE/SALT2, and the code developed by iHoweU et all 
( l2005i) . Difficult spectra were discussed, on a case-by-case basis, 
until agreement was reached. In case of disagreement, the most 
conservative typing was chosen. This procedure makes the iden- 
tification of the SNe la presented in th is paper homogeneo us to 
the identificat i on of Gemini spectra of iHowell et al.l (l2005h and 
iBronder et al.l ( 1200 8'). Note however that we do not use here the 
confidence index (CI) classification of Howell et al. (2005). As 
a guide, the SNe la of the present paper correspond to CI=4 and 
5 while our SNe la* correspond to CI=3. 

5. Results 

5. 1 . Individual spectra 

In this section, we present the spectra of the 124 identified 
SNe la of the SNLS 3rd year VLT spectroscopic survey, along 
with their identification as SN la, SN la* or SN Ia_pec. Only 
two objects have been identified as SN Ia_pec (SNLS-03D4ag 
and SNLS-05Dlhk). Their properties are described individually 
in Section |53] In the following, we include them into the SN la 
subs ample. 

Table[T]lists the SN la and SN la* spectra and their observ- 
ing conditions. In some cases, several spectra of the same can- 
didate were taken due to poor conditions at the telescope, or an 
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insufficient S/N for a secure identification. The asterislc in the 
coiumn Dsep denotes cases for which a separate extraction of the 
SN from the host component was not possibie (inciuding, for 
some candidates, the absence of a detectabie host). 

Tabic |2] summarises the resuits of our identification for each 
of the 139 spectra corresponding to the whole set of 124 SNe la. 
Name, type, redshift, redshift source, i.e. host (H) or SN (S), 
phase, host type, fraction of host used in the best-fitting model, 
and average S/N per 5 A bin are given. This latter quantity is 
computed on the host subtracted spectrum (when relevant) for 
each SN. PEGASE best-fitting host model is given as a letter 
for the morphology (E, SO, Sa, Sb, Sc, Sd) followed by a fig- 
ure between parenthese s indicating the age ( in Gyrs). When the 
best-fit is obtained for a lKinnev et alj (1 19961) template, we indi- 
cate the two contiguous Hubble types between which the best- 
fitting galaxy model is interpolated. Note that a host component 
is always allowed (this includes a null template), even when the 
SN spectrum has been extracted separately from its host. Where 
the best-fit is obtained for a model with no host contribution, the 
label 'NoGalaxy' is given. 

Table [3] gives the number of SNe for each class (SN la and 
SN la*) per field and in total. Figures [ATT] to |A. 1 3"9l show the 
full (host+SN) PHASE extracted spectrum (left panel). For each 
spectrum, the name of the SN is followed by the date of spec- 
troscopy, in terms of the number of days elapsed since January 
1st, 2003. The red dashed line is the SALT2 model obtained with 
no recalibration, the red solid line being the same model after 
recalibration. Whenever fitting with a host galaxy template is 
necessary, we show the best-fitting host spectrum model (blue 
solid line). In these cases, we present in the right panel the host 
subtracted SN spectrum obtained by subtracting the host model 
(blue solid line in left panel) from the PHASE extracted spec- 
trum. The best-fit, recalibrated, SN la model (solid red line) is 
overplotted on the host-subtracted SN spectrum. 

As discussed above, host subtraction is a key issue and tech- 
niques are not perfect. In particular, the PEGASE templates we 
use have no emission lines. As a consequence a number of our 
spectra show residual host lines, e.g., SNLS-04Dlsa (Fig. |A.40] i 
with residual [O ii] emission and Ca ii H&K absorption. In some 
cases, the fit is poor in some portion of the wavelength range. 
There can be different explanations for this. In the case of SNLS- 
04Dlhx (Fig. IA.27I ). there are two galaxies along the line of 
sight. The PHASE host model used for the extraction is not accu- 
rate, and the extracted spectrum shows strong host residuals. In 
the ca se of SNLS -03 D4gf, SNLS-03D4gg and SNLS-04D2cw 
(Figs. |A191Ia;201 and |A32] i, the SALT2 model is unable to re- 
produce the UV wavelength region due to the lack of UV cover- 
age in the SALT2 training sample used. More specific comments 
on individual SNe la are given in the corresponding caption. 

5.2. Average properties of the SN la and SN la* samples 

The main parameters characterizing the SN la and SN la* sub- 
samples are given in Table H] and are discussed below. Figures [T] 
and |2] show the redshift and phase distributions of both SN la 
and SN la* samples. Redshifts range from 0.149 to 1.031. 
As expected, the SN la* subsample has a higher average red- 
shift (< z >SN la*- 0.70 ± 0.03) than the SN la subsample 
(zsN la - 0.60 + 0.02). The average redshift of the whole sample 
is < z >= 0.63 and the median redshift is 0.62. Below z ~ 0.3, 
all SNe la are identified as certain SN la. 

The average phase of the SN la* subsample is significantly 
higher (< >sn la*- 5.8 + 1.1 days) than for the SN la 



(< (f) >SN la- 1.1 +0.4 days). This is generally caused by similar- 
ities between SN la spectra one week past maximum and SN Ic 
spectra. These cases are labeled SN la* unless Sin is clearly 
seen. Fig.|2]shows that most spectra at phases later than 10 days 
are classified as SN la*. This also reflects the lower-S/N in these 
(fainter) spectra. The average S/N per 5 A bin is ~ 6 for SN la, 
~ 2 for SN la*. 

We now compare the rest frame B-band magnitudes, at max- 
imum light, for the two samples. For each SN, m^, the ap- 
parent rest frame B-band magnitude, is determined as part of 
the light curve fit with SALT2. In order to assess the signifi- 
cance of any discrepancy between the SN la and SN la* sam- 
ples in terms of physical properties of the SNe, we compute, 
for each SN, the "distance corrected" magnitud^l m*'" = m* - 
5log(c^^ Hodiiz,®)), where is the luminosity distance, c is 
here the speed of light and = {Ho, Q^, ^a] the set of cosmo- 
logical parameters of the underlying c osmology. Here , we adopt 
the values = {70,0.27,0.731] (.Astier et alj|2006h . We find 

< >sNia= 23.95+0.04and< m*^ >sn la*^ 24.12+0.08. As 
expected, the SNe we classify as SN la* are fainter, on average, 
than the ones we classify as SN la (by < Am*'" >= 0.17 + 0.09 
mag). This reflects the fact that low S/N spectra at a given phase 
fall preferentially into the SN la* category, as they are more dif- 
ficult to identify. 

It is of interest to account for the difference found in the 

< m*g > values. One obvious explanati on is that it i s re- 
lated to the well -known "brighter- slow er" (|Phillipslll993l) and 
"brighter-bluer" ("Tripp & Branch' '1999*) relationships observed 
in SNe la populations. Figures [3] and |4] show the SALT2 colour 
c and xi parameter distributions for the two subsamples. The 
SN la* appear to be slightly redder, as their average colour 

< C >SN la* — 0.05 ± 0.03 is higher than for SN la: < c >sn la- 
0.009 ±0.013. Regarding xi, we find <jci >sn la^ -0.02 + 0.09 
and < xi >sN la*- 0.004 + 0.18, which are consistent within 
Icr errors (the SN la sample having faster, narrower light curves 
than the SN la* sample, on average). The "brighter-slower" and 
"brighter-bluer" relationships translate into a magnitude differ- 
ence between the two samples of < Am^ > - -a^ < Axi > 
+/3 < Ac >. Using a, = 0.13 and /? = 1.77 jGuv et al.ll2007h . 
we find < Am«' >^ 0.07 ± 0.06. The "observed" < > dif- 
ference between the two samples is therefore consistent with the 
"brighter-slower" and "brighter-bluer" relationships. This gives 
us confidence that our SN la* sample is not significantly pol- 
luted by the inclusion of red, non-SN la objects, such as SNe Ic. 
We can estimate the contamination of the SN la* sample by 
core-collapse SNe by comparing the number of certain SN Ib/c 
we have typed to the number of certain SNe la. In the whole 
SNLS data, ~ 3% of securely identified SNe are of the SN Ib/c 
type. Applying this ratio to the 38 SN la* of the VLT sample 
shows that we should expect at most 1 contaminant. Note how- 
ever that the detection efficiency of the SNLS does not allow us 
to detect SN Ib/c SNe beyond z < 0.4 (Bazin et al. 2009). As we 
only have 3 SN la* with z < 0.4, this leads to a more realistic 
estimate of 0. 1 SN la* that might be contaminant. 

5.3. Peculiar SN la 

We have identified two SN la in our sa mple showing strong simi- 
larities with the spectra of SN 1999aa jGaravini et al.ll2004l) (see 
Fig. EJ. These are 03D4ag and 05Dlhk (this latter SN was al- 

' Note that m'^' is not the quantity used to constrain the cosmological 
parameters, as it is not corrected for xj and c. 
* Ho is in units of km/s/Mpc. 
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ready noted as a peculiar SN la by lEllis et aT]|2008l) . Table |5] 
summarises the properties of these two SNe as well as the param- 
eters obtained when fitting their spectra with SALT2. Both SNe 
are at relatively low redshift (z = 0.285 and 0.263 for 03D4ag 
and 05Dlhk respectively) and at early phase (~ -9 and ~ -5 
days). 05D 1 hk has a large xi value, which translates into a Amis 
parameter ("Phillips" 1993) of 0.89 using the formula given in 
[(juy et al (2007)^ By comparison, SN 1999aa has Amis = 0.85 
Ijha et al.ll2006h . For 03D4ag, we find Amis = 0.97. Visual in- 
spection of the spectra shows that 03D4ag and 05Dlhk have 
properties typical of high stretch SN la: shallow silicon absorp- 
tion and a blue spectrum, even for the early phases under con- 
sideration. 

The possibility of SNLS-03Dlco being a 1991bg-like event 
has been discussed in iBronder et al. I (12008 '). as there is also a 
Gemini spectrum of this event. iBronder et al. (2008) found that 
03Dlco had a large Mgii A 4300 equivalent width (EW), con- 
sistent with the excess absorption measured in under-luminous, 
low-z SNe spectra. However, the presence of Si ii A 4000 in the 
Gemini spectrum contradicts this hypothesis, as this feature is 
not seen in under-luminous SNe because of extra absorption due 
to Tin and Feii in the same wav elength range. Based on this 
observation. iBronder et al.l (l2008h concluded that 03Dlco was 
likely t o be a "normal" SN la, in agreement with its normal light 
curves (" Astier et al.ll2006l) . The SALT2 fit of the VLT spectrum 
of 03Dlco confirms this conclusion (see Fig. |A.5t . 

5.4. Contamination by non-SN la 

An obvious challenge in the identification of distant SNe la is 
to avoid confusion with other types, such as SNe Ic, especially 
one week after maximum light and beyond. Template cross- 
correl ation techniques such as SNID can help dMatheson et al.l 
l2005b . but still rely on (unavoidably) incomplete template li- 
braries. Our SALT2 identification procedure provides some ad- 
ditional information on non SN la types (or peculiar SNe la) 
through the fitted parameter values, which can differ from those 
of a typical SN la, e.g., an unusually red colour, and/or a high 
xi parameter, and/or a high recalibration parameter This is not 
equivalent to a direct identification of non SN la objects, but does 
allow a mechanism by which peculiar events can be identified. 
We explore these parameters in this section. 

We start with the SALT2 "recalibration parameters" used to 
adjust the SALT2 photometric model to the observed data. We 
focus on the first order recalibration parameter yi (Section |4~2] l. 
This parameter can be interpreted as the "tilt" required to ad- 
just the observed spectrum to the SALT2 colour model. A large 
recalibration is a sign that the SALT2 model is not able to prop- 
erly model the data, as would be the case if we were trying to fit 
a non-SN la spectrum, but could also occur for a SN la spectrum 
whose properties were very different from those of the training 
sample. Large tilts can also be needed when the host subtraction 
has failed due to inadequate host modelling or too strong a con- 
tamination, although no strong correlation exists between yi and 
the host fraction. 

The average yi value for the SN la* is < yi >= 0.26 ± 0.29, 
with < yi >= 0.26 + 0.11 for the SN la: the required tilts needed 
to recalibrate the photometric model are only moderate. The dis- 
persions in the yi values are quite large for both samples, with a 
larger variation for the SN la*: cr^^^ = 1.1 and cr^f^ = 1.9. 
This can be explained by the inclusion in the SN la* sample 
of more host-contaminated spectra. If we select only SNe for 
which a separate extraction from their host was not possible, the 
mean host fractions fg (i.e. the contribution of the host model 



to the full spectrum averaged over the whole spectral range) are 
< fg >= 0.32 + 0.03 and < fg >= 0.47 + 0.05 for the SN la and 
SN la* subsamples, respectively (see Table |4|i. Clearly, spec- 
tra identified as SN la* are more host contaminated than SN la 
spectra. 

We do not find any systematic trend with redshift or phase 
in the "tilt" values needed to accomodate the SALT2 model with 
the spectra. In Fig. |6] we show the magnitude Am of the recal- 
ibration of the photometric model as a function of wavelength, 
for a subset of our spectra chosen between z - 0.4 and z = 0.6 
(< z >^'v/fl^ 0.512 and < z >s'^^«*= 0.518 for this subset). 
The light blue area is for SNe la*, the dark blue one SNe la. 
At this redshift, we find that a « 15% recalibration is needed at 
both ends of the effective spectral range, for both the SN la and 
SN la* categories. 

A potentially more quantitative criterion is the (reduced) 
xl of the SALT2 spectral fits. We show in Fig. [T] the reduced 
xl histogram of all SALT2 fits (both for the SNe la and 
SNe la* samples). The bulk of the xl values are around 1 
(or slightly higher), but with a tail of 1 1 objects with xl > 
(Fig. Ill the two highest xl objects are not shown). With these 
objects, we find < xl >- 1-33, while excluding them yields 
<xl >- 1.11. Among the 11 high;^^^ objects, one finds 03D4ag 
and 05Dlhk, the over-luminous events identified in Section [53] 
Their best-fit xl values are 2.48 and 3.12 respectively. 
The other objects are (ordered by increasing xl)- SNLS- 
03Dlfc(2.01), SNLS-04D4ht(2.03), SNLS-05D2ct(2.04), 
SNLS-05D4cw(2.18), SNLS-04D2bt(2.23), SNLS- 
04D2fs(2.42), SNLS-04Dldc(3.24), SNLS-03D4au(4.63) 
and SNLS-05D4fr(16.33). 

03Dlfc (Fig.lAJl SN la at z = 0.332) has a separate extrac- 
tion from its host. The spectrum has a high S/N and the SALT2 
fit parameters are typical of a normal SN la. Si ii /i 6150 is visible 
but shallow. The spectrum is not as blue as the one of 05Dlhk 
(which is at about the same phase, i.e. -5 days), but the photomet- 
ric model requires three recalibration parameters to consistently 
fit the spectrum. Although the stretch is normal (5= 1 .006), this 
SN might be slightly peculiar 

04D4ht (Fig. rATTl SN la at z = 0.217) is heavily host con- 
taminated (> 60% of host in the best-fitting model). It has an 
high colour value (c - 0.5). 04D4ht is very red but with strong 
Si II in its spectrum and is identified as a SN la. However, it is 
located very near the core of its host (a late type spiral) and po- 
tentially heavily extinguished. 

05D2ct (Fig. |A!95l SN la* at z = 0.734) is another example 
of a SN close to its host centre (more than half of the extracted 
signal is modeled by a Sd galaxy), slightly red (c - 0.14) but 
with a normal stretch (s = 0.994). It is likely to be extinguished. 
Due to its fairly high redshift and host contamination, no Si ii is 
visible either at 6150 or at 4000 A and it is identified as a SN la*. 

05D4CW (Fig. IA.116I SN la at z = 0.375) is a blue (c = 
-0.15) SN. It is heavily host-contaminated, probably an early 
type galaxy. As it is close to its host centre, a separate extraction 
was not possible. However, Sin A 6150 is clearly visible, and it 
is classified as a SN la. 

04D2bt (Fig. [A46l SN la at z = 0.220) is located in the 
bulge of its early type host. It is red (c = 0.18), likely due to 
host extinction. Si ii and S ii are clearly visible, the classification 
is SN la. 

04D2fs (Fig. |A34l SN la at z = 0.357) has a large S/N with 
Si II and S ii clearly visible. The high xl value can be explained 
by the high S/N (low noise) of the spectrum. 
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04Dldc (Fig.|A25j SN la at z = 0.21 1) is similar to 04D2fs 
with a very high S/N. Si ii and S ii are obvious. Again, the high 
xl value reflects the small errors on the spectrum flux. 

03D4au (Fig. |Al2l SN la* at z = 0.468) is red (c = 0.18) 
and host contaminated, a week past maximum light. It is located 
right at the centre of its late-type host, probably extinguished, 
and with strong emission lines that are difficult to subtract, ex- 
plaining the high;^'^ value. 

05D4ff^ (Fig. IA.123I SN la* at z = 0.402) is heavily buried 
in its late-type host. The presence of strong [O ii], [O iii] and H/3 
emission lines, not present in the host PEGASE template, ex- 
plains the poor;i'^. 

In conclusion, a cut in xl can help with the identification of 
pecuUar (and possibly non-SN la) events, but it is not a straight- 
forward identification. Several effects can conspire to give a high 
xl value, even in the case of obvious SNe la. The results of this 
section suggest that our two subsamples are unlikely to be sig- 
nificantly contaminated by non-SNe la. 

6. Composite spectra of the SN la sample 

In this section, we build average VLT spectra in six regions of 
phase-redshift space to compare the evolution of their properties 
with redshift. 

6.1. Methodology 

We select spectra with phase <p < -4 days (pre-maximum spec- 
tra); -4 < < 4 days (maximum spectra); (p > 4 days (post- 
maximum spectra) for SNe at redshift z < 0.5 and z > 0.5, 
respectively. We only use SN la spectra, as the SN la* spec- 
tra number in each of the six regions is small and even zero 
for z < 0.5 around maximum (since the spectra in this phase- 
space region are of sufficient S/N and quality to be identified as 
SN la). The number of spectra in each phase bin for both z < 0.5 
and z > 0.5 is summarised in Table |6] For the SN la spectra, 
we have 12 pre-maximum spectra (of which 7 are at z < 0.5), 
57 spectra at maximum light (of which 12 are at z < 0.5) and 
24 post-maximum spectra (of which 6 are at z < 0.5). The av- 
erage redshift is < z >= 0.36 + 0.02 for the z < 0.5 sample 
and < z >= 0.70 ± 0.02 for the z > 0.5 sample. We have more 
spectra in the z > 0.5 bins (68 spectra) than in the z < 0.5 bins 
(25), as the median redshift of our sample is 0.62. The spec- 
tra of 03D4ag and 05Dlhk are excluded from the z < 0.5 pre- 
maximum light region, as they have been identified as peculiar 
SNe la (see Section lSJl i. 

To construct the average spectrum in each region, all indi- 
vidual spectra (shown in the right panels of Fig. IA.ll to |A.139l 
whenever a host model has been subtracted) are brought into the 
rest frame and rebinned to 5 A. All spectra are colour-corrected 
using an up-to-date version of the SALT2 colour law and the 
colour value c obtain ed from the SALT2 fit for each SN. Using a 
ICardelh et all (1 19891) extinction law with = 3.1 and E{B - V) 
values obtained for each SN from their light curve fit, yields 
very similar results in the wavelength range under consideration. 
We adopt the SALT2 colour correction in the following. Fluxes 
are normalised to the same integral in the range 4450-4550 A. 
For each wavelength bin, an average weighted flux and its cor- 
responding uncertainty are computed from all spectra included 
in this bin. Because of the variety of the redshifts involved, the 
number of spectra entering the average varies from one bin to 
the other, as does the average phase in a given bin. In prac- 
tice, the number of spectra in each bin decreases at both ends of 



the wavelength scale but is approximately constant in the range 
4000-7000 A. 

As a consequence of this averaging procedure, the mean 
phase of the average spectra varies from one wavelength bin to 
the other. In practice, this variation distorts only very moderately 
the spectra as the mean phase is roughly the same from one end 
to the other of the wavelength scale. However, for a given phase 
range (pre-maximum, maximum or post-maximum), the mean 
phases of the average spectrum at z < 0.5 can significantly differ 
from the corresponding one at z > 0.5, which should be borne in 
mind when comparing the average spectrum at z < 0.5 with that 
at z > 0.5. In practice, due to a similar phase sampling of our 
spectra with redshift, this difference is marginal. We find the fol- 
lowing average phases (in days): (/>^^"'^ = -6.1 and<p'^^-^ = -5.4 
for pre-maximum spectra, 0^*"-'' - 0.5 and (p^^*^-^ - 0.5 for spec- 
tra at maximum, and 0^^"-^ = 7.7 and <p~^^-^ =5.5 for post- 
maximum spectra. 

6.2. Comparison to Hsiao et al. template spectra 

Figure[8]shows the result of averaging the SN la spectra (in blue). 
The left column is for z < 0.5 spectra, the right for z > 0.5 
spectra. From top to bottom, pre-maximum, maximum and post- 
maximum s pectra are shown. Average spectra built from the 
iHsiao et aTl d2007 ) template series (red curve) are overlaid on 
top of the VLT average spectra. N ote that about two -thirds of 
the spectra used in constructing the iHsiao et al.l (l2007h template 
are at low redshift (z < 0.1) and lack UV coverage, so the UV 
section of the red curves come from the remaining high-redshift 
spectra used in building the Hsiao et al, (.2007) template. In each 
region of Fig. |8] phases in each wavelength bin are the same 
for the template and the VLT composite spectrum, weighted in 
exactly the same way. 

The overall a greern ent between the VLT average spectra and 
the iHsiao et aTl (l2007l) template is good in all region s. Note 
that n o colour correction has been applied to the Hsia o et all 
(|2007|) templates. As this template is designed for use in light 
curve fitters that im plement "warping" techniques (e.g., SiFTO, 
IConlev et al.ll2008h . there is no a priori reason that the contin- 
uum should agree with our composite spectra. However, we find 
that the agreement is almost perfect in the optical region of the 
spectra except around the Ca ii A 3700, Si ii A 4000 and Si ii A 
6150 features, for the maximum average spectrum at z < 0.5. 
In the UV region, one notices more discrepancies in the fluxes 
(a low number of spectra are used in this region). Part of this 
effect might be due to the fact that we normalise the spectra in 
the optical region. This spectral region is also the most sensi- 
tive to differential slit losses. We find a satisfying match of the 
positions of the UV A 3000 - 3400 features. Note that the peak 
around 3200 A decreases from pre- to post-maximum phases in 
a proportion that is well reproduced by the iHsiao et al.i (i2007l) 
model, both at z < 0.5 and z > 0.5. 

6.3. Comparison ofz < 0.5 and z > 0.5 spectra 

We now compare the average spectra at z < 0.5 and z > 0.5 for 
pre-maximum phases (upper panels of Fig. |9]l, maximum-light 
(middle panels) and post-maximum (lower panels). The compar- 
ison is done in the region of intersection, from the UV up to the 
mid-optical wavelengths. For each panel, the blue curve is the 
z > 0.5 average spectrum, and the black curve is the correspond- 
ing z < 0.5 spectrum. At the bottom of each panel, we plot the 
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residual A^^'^(A) = f-^\A) - f^H^) (solid black line) with 
+ lcr errors. 

The mean spectra can be different for two reasons: they 
might be intrinsically different (i.e. evolution in spectral prop- 
erties), or the phase distribution of the samples can differ. We 
theref ore compute the res idual A^*""'(/l) (solid thick red curve) 
for the lHsiao et alj (1200 7') templates shown in Fig. [8] As the un- 
derlying assumption in building such spectral templates is that 
there is no evolution between low and high redshift, A^^'"°(A) 
should in principle be zero over the whole spectral range. Any 
deviation from zero, in a given wavelength range, can be at- 
tributed to a difference in the average phase of the z < 0.5 and 
z > 0.5 composite spectra in this range. When inspecting differ- 
ences in the VLT composite spectra at z < 0.5 and z > 0.5, it 
is important to refer to A""'"'(^): if A^^^(^) follows A"'''"'(A), 
deviations can be traced back to differing phase distributions of 
the z < 0.5 and z > 0.5 composite spectra. In the opposite case 
^^vir jjQj. 2gj-Q while A"'""°{A) is zero), any differences should 
be real. 

To quantitatively assess the significance of the differences, 
we define a reduced measure of the agreement of the z > 0.5 
VLT composite spectrum and the z < 0.5 composite spectrum. 
The variance entering the definition of xl is the sum of the vari- 
ances of the two composite spectra. A value ~ 1 indicates that 
the two spectra are consistent with one another We now examine 
each phase bin in turn. 



iBronder et alj (l2008l) indicate a possible difference in the 
EW of Mgn A 4300 between low and high redshift, but they 
conclude that it is likely due to differences in the epoch sam- 
pling and number of obje cts at low- and hig h-redshift. This dif- 
ference is not found by iFolev et akl (l2008a) when comparing 
their composite ESSENCE spectrum with a Lick low-z com- 
posite (note that they use a simpler method for host subtraction, 
which may alter the signific ance of their comparison). More re- 
cently, [SuU^aneTal] ( l2009l) studied the possible evolution in the 
EW of intermediate mass elements (IMEs) in the < z < 1.3 
redshift range. They do not have a Mgii measurement in their 
highest-redshift bin, but the predicted variation of Mg ii EW is 
consistent with zero over their redshift range. 

The Si ii-Fe ii-Fe iii A 4800 blend feature is shallower in our 
z > 0.5 s pectrum than at z < 0.5, in qualitative agreement with 
iFolev et al. (2008a), who find that this feature is much weaker in 
their ESSENCE spectrum. They interpret this difference as due 
to a weaker Fe iii /I 5129 line in the high-z spectrum. 

The Sill A 4000 feature is shallower at high re d shift, in agree- 
ment with the findings of ' Sullivan et al.l {200^. 'Bronder et alj 
(2008) do not mention such a difference, but Foley et al. (2008^ 
find the same kind of trend for Sin A 6150 (this feature be- 
ing shallower in their high-redshift ESSENCE spectrum than in 
the low-redshift Lick spectrum), although they do not mention 
this for Si n A 4000. We do not find this difference in our pre- 
maximum spectra, though it may be present at post-maximum. 



6.3.1. Pre-maximum spectra 

For pre-maximum spectra (0 < -4), the composite spectra are 
typically consistent (xl ~ 1.14 for 588 degrees of freedom - 
D.O.F). The spectra are most discrepant in the UV region (3300 
to 4000 A). A possible explanation is that, as the bands bluer 
than the rest frame B-band are given a lower weight in the 
SALT2 light curve fit, the colour correction applied to the spec- 
tra that uses the colour parameter c derived from the light curve 
fit, is more efficient in the red than in the blue. Alternatively, 
this might reflect a g reater variability in this spectral region 
(e.g.. lEllis et al.ll2008l) . More variability is expected here due to 
variations in metal licities of the progen i tors, particularly at pre- 
maximum phases dHoeflich et all [19981: lLeinzetan '2()Q0). The 
residuals in Fig.|9]show the UV region of the pre- and maximum 
spectra shows more discrepancy than at post-maximum. 

6.3.2. Maximum-light spectra 

The largest discrepancy between z < 0.5 and z > 0.5 is found 
at maximum light (-4 < <p < 4). This is the phase where we 
have the highest number of spectra (12 at z < 0.5, 45 at z > 0.5), 
and the statistical errors are the smallest. We find;^'^ ^ 4.80 for 
601 D.O.F. - the two average spectra are formally not consistent. 
The largest differences are seen around Ca n A 3700, Si ii A 4000, 
Mg II A 4300 and Fe ii A 4800. Overall, the z < 0.5 maximum 
spectrum has deeper absorptions than its z > 0.5 counterpart. 

In the UV up to 3500 A, the residuals A^^'^(A) correlate with 
the template residuals A^-^"'°(A), and the discrepancies between 
z < 0.5 and z > 0.5 may be due to differences in mean phases 
between z < 0.5 and z > 0.5 spectra (recall that due to our 
averaging procedure, the mean phase varies with A and depends 
on the phase of the spectra used to build the average spectrum 
in a given wavelength bin). Around Ca ii A 3700 and beyond, the 
discrepancies are seen in the spectra, not in the template, and 
could be real. 



6.3.3. Post-maximum spectra 

At post-maximum {(p > 4), we find a good agreement be- 
tween the z > 0.5 average spectrum and its lower z counter- 
part: xl ~ 1 10 (547 D.O.F.). Variations in the spectral residuals 
closely match those in the template residual and reflect mean 
phase variations between the z < 0.5 and z > 0.5 average spec- 
tra. Once again, the level of discrepancy, though small, is highest 
in the UV and in the Si ii A 4000 region. It is interesting to note 
that, as for the maximum light spectra, the flux at the position of 
this line (which has almost disappeared a week past maximum) 
is shallower for z > 0.5 than for the z < 0.5 spectrum. For Mg ii 
A 4300 and Fe ii A 4800, no difference is found with redshift. 

Inspection of the spectra shows the presence of a residual 
[On] host line in the z > 0.5 spectrum, for which a separate 
extraction of the SN and the host is often difficult. Moreover, 
our PEGASE host templates do not have emission lines. Other 
residual absorption (e.g. Can H and K and Balmer lines) is also 
found in individual spectra. 

Our main conclusion of this redshift comparison is that de- 
spite the overall agreement between low (z < 0.5) and high 
(z > 0.5) spectral data, some discrepancies are found in char- 
acteristic absorption features. This may indicate evolution with 
redshift, or a signature of some selection effect. We discuss this 
in the next section. 



7. Discussion 

In recent years, various large-scale SN programs have published 
sets of SN la spectra at intermediate (Balland et al. 2006, 200^ 
Zheng et al. 2008; Foley et al. 2008b) and distant (Howe Uet al 
2005; Lidman et al., .2005; .Matheson et al.1 12008^ .Bronder et af. 
2008; Ellis et aLl l2008HFolev et al.ll2008al) redshifts. The SNLS 



VLT SN la spectral data set presented in this paper supplements 
these existing sets with 124 new SNe la, or probable SNe la. 
This constitutes the largest high-redshift SNe la spectroscopic 
sample published so far 
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Comparing our d ata to this lit e rature , we have the highest 
SN Ia*/SN la ratio. ' Howell et all (l2005h publish the 1st year 
of SNLS spectroscopy at Gemini. The redshift range targeted 
is 0.155 < z < 1.01, the median redshift being 0.81, higher than 
our median redshift 0.62. Attributing an index of spectral qual- 
ity to the Gemini spectra, they identi fy 34 SN la and 7 S N la* 
among 64 SN candidates (recall that Howell et al.l (|2005[) index 
CI =4 and 5 to SN la and index Cl=3 to SN Ia*). lLidman et alj 
( l2005h find 15 SN la and 5 SN la* in the range 0.212 < z < 
1.181 (< z >«; 0.5). In both cases, SN la* amount to -20% of 
their total SN la sample. We find 86 SN la (of which two are 
SN la_pec) and 38 SN la* in our VLT s ample, that is ~ 30 % of 
SN la*. Ot her surveys such as SDSS-II (IZhe ng et alj|2008l) and 
ESSENCE dMatheson et al.l ! 2005) both find about 10% of prob- 
able SN la (equivalent to our SN la*). The SDSS-II survey tar- 
gets low to intermediate redshifts (0.05 < z < 0.4, < z >~ 0.22), 
which probably makes identification easier The average redshif t 
of the ESSENCE SNe la is < z >~ 0.4 jFolev et alj|2008bl) . 
significantly lower than ours. 

We identify fewer SN la* towards the end of the SNLS 3rd 
year period. The average spectroscopic MJD of SN la classified 
spectra is 53393 whereas it is 53272 for SN la* spectra. As all 
spectra are treated on an equal footing from the point of view 
of extraction and identification, this can not be due to improve- 
ment in data processing or refinements in classification, but more 
likely indicates an improved target selection efficiency and opti- 
misation of spectroscopic time during the course of the survey. 

A key result of this work is the construction of average spec- 
tra for various phases and redshift bins from the homogeneous 
sample of VLT SNe la. We have used these high quality average 
spectra to compare spectral properties at z < 0.5 and z > 0.5. 
We find differences in the depth of some optical absorption fea- 
tures (Sin and Can) around maximum light. We find that the 
absorptions due to these elements in the z < 0.5 spectrum have 
weakened in t he co rresponding z > 0.5 spectrum. Recently, 
ISullivan et al.l (l2009t) have shown evidence for a similar weak- 
ening of singly ionize d IME EWs a t highe r redshift. Using data 
from various searches. 'How ell et al. |(|2007) find an ~8% increase 
in average light curve width for non-subluminous SNe la in the 
0.03 < z < 1.12 redshift range and interpret it as a demographic 
evolution. High-z SNe la have an average stretch higher than 
at low redshift, thus are more luminous and hotter, and should 
ionize more IMEs, depleting singly ionized absorptions in high 
redshift spectra (e.g., Ellis et al. 2008). 

The average stretches of our two redshift subsamples are 

< s >'<°^= 0.975 + 0.017 and < s >'^0-5= 0.983 + 0.009 re- 
spectively. These two values are similar gi ven the uncer t ainties . 
This is not inconsistent with the results o f iHowell et al] (|2005|) : 
ISullivan et al.i (2009), as the predicted change in the spectra over 
our redshift range due to a demographic shift is likely to be very 
small. The average colour^ are < c >'<°-^^ 0.06 + 0.03 and 

< c >^-°-^= -0.01 + 0.01. The average rest frame distance cor- 
rected magnitude s (using = {//o,Qm, ^^aI - {70,0.27,0.73} 
lAstieretal.ll2006l) ai-e < m*" >,-<o.5= 24.10 + 0.10 (cr = 0.52) 
and < ml' >z>o.5= 23.88 + 0.03 (cr = 0.26), the z > 0.5 
subsample having brighter supernovae, on average, in qualita- 
tive agreement with what is expected from Malmquist bias, as 
brighter supernovae, preferentially selected at higher redshift, 
tend to be bluer. In order to assess the significance of the ap- 
parent spectral differences in our maximum light spectra, we 
have built two new subsamples (one for each redshift bin) by 



' We limit our set to SNe la with c < 0.3 to reject red, heavily extin- 
guished objects. This cuts three SNe la out of the sample. 



selecting spectra with phase in the range -4d < (p < Ad and 
bluer than average (i.e. spectra of SN la with colour c < 0). 
Using these cuts, we end up with only 6 spectra for the z < 0.5 
bin, and 26 spectra for the z > 0.5 bin. As we expect to ob- 
serve bluer objects at higher redshift due to Malquist bias and 
to the 'brighter-bluer' correlation, this should select two sub- 
samples with roughly comparable photometric properties. The 
average distance corrected magnitudes and dispersions of these 
new subsamples are < m*^ >;<o.5= 23.77 + 0.04 (cr - 0.11) 
and < m*^ >z>o.5= 23.75 + 0.03 (cr = 0.17). The two first mo- 
ments of the m*g distributions are thus similar. We build new 
average spectra at maximum light in the same way as in Section 
|6] The results are shown in Fig. [10] and must be compared to 
the middle panel of Fig. |9l Most of the differences have diseap- 
peared, except around 3700 A and 3900 A. We find xl = 0.97 
(544 D.O.F). The difference at 3700 A is due to a residual [O ii] 
line in the z > 0.5. If the differences seen in Fig. |9] were due 
to imperfect calibration, slit losses or imperfect host subtraction, 
we would expect to see comparable differences when selecting 
the two 'blue' populations shown in Fig. [10] The fact that the 
composite spectra nicely overlap for these two populations is in 
itself an indication that the differences with redshift in our maxi- 
mum light spectra are real. We have checked that the differences 
remain when selecting a subset of spectra with a host galaxy 
fraction below 20%. They also remain when we select a subset 
of spectra with little recalibration (lyil <0.5). Thus it is unlikely 
that the differences are due to slit losses or a poorer host subtrac- 
tion at higher redshift. These differences more likely result from 
the selection of brighter and bluer supernovae at higher redshift. 

We now com pare our average spectra to the ones of 
Ellis_etalJ (|2008|) . These were built from spectra obtained at 
Keck as part of detailed studies programs using SNLS SNe. 
Some SNe la ar e common to both samples. We find a qualitative 
agreement with lEUis et al.l (l2008b regarding variations in the UV 
part of the spectra. Figure [TTj shows the comparison between a 
subset of 51 VLT average spectra between z = 0.35 and z = 0.7 
(with < z >~ 0.5 for this subset), for pre- and maximum phases 
(in black) with the average spectra obtained from 23 high S/N 
Keck spectra, with < z >~ 0.5 (in red, from Tables 3 and 4 of 
their paper). Both VLT and Keck average spectra are shown with 
lower and upper 2cr errors and have been normalised to unity in 
the same wavelength region (at 3900 A). VLT spectra have been 
recalibrated using the SALT2 recalibration function at first order 
(the "tilt" coefficient y\). The average phases of the Keck spectra 
are comparable to the VLT ones. 

The two sets of pre-maximum spectra are consistent: xl ~ 
1.00 (500 D.O.F). For around maximum spectra, we find;^^^ 
1.17 (500 D.O.F). The VLT average spectra are on the lower 
limit of the Keck spectra in the UV region for both phase ranges. 
In the optical region, the average spectra are in remarkable 
agreement for pre-maximum spectra. For maximum-light spec- 
tra, the agreement is fair in the optical, with some structure seen 
in the residuals. It is in the UV region that differences are most 
noticeable. There could be several explainations: this region is 
more sensitive to differential slit losses, and we have normalised 
the spectra at 3900 A (had we chosen a shorter wavelength, the 
UV discrepancies would have been less apparent, but the optical 
differences greater). Also, the SALT2 col our law has been used 
to build the average VLT spectra, whereas [Ellis et al" I (120081) use 
the SALTl colour law, though this has a negligible effect on the 
average spectra. We find that the Ca n /I 3700 and the Si ii A 4000 
absorptions are consistent, both in pre-maximum and at maxi- 
mum spectra. 
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Recently, iFolev et aTl (l2008al) have published composite 
spectra from ESSENCE data and they compare them in various 
bins of redshift and phase, to the Lick z < 0.1 counterpart com- 
posite spectrum. Comparing the ESSENCE and VLT spectra in a 
similar fashion to our comparison with the Ellis spectra, we find 
sub stantial differences a t all phases. This may be due to the fact 
that lFoley et alJ (l2008ah subtract a single average host spectrum, 
derived from a PCA analysis, from all of their SNe. Regarding 
specific spectral features, they find that the most important dif- 
ferences at maximum light between low- and high-redshift are 1) 
Fe 11 ^ 5 129, that yields a smaller EW of the {Fe ii A 4800} blend 
in the ESSENCE spectra than in the Lick composite, and 2) the 
lack of absorption at 3000 A in the highest redshift bins of their 
ESSENCE composite spectra. We qualitatively agree with their 
analysis of Fe ii A 4800. We however find the presence of absorp- 
tion at 3000 A in our average spectra, for all phases at z > 0.5, 
as seen in Fig. |8] For z < 0.5, the absorption is at the limit of our 
effective spectral range (see Fig. [8]), but the absorption seems 
present, at least for pre-maximum and at maximu m spectra (it 
is very noisy for post-maximum spectra). Following iFolev et alj 
(l2008ah . this difference could arise from a selection bias in the 
ESSENCE sample that is not present in the VLT sample. 



8. Conclusion 

We have presented 139 spectra of 124 SNe la (SN la type) or 
probable SNe la (SN la* type) at 0.149 < z < L031 ob- 
served with the FORSl instrument at the VLT during the first 
three years of the SNLS survey. This is the largest SN la spec- 
tral dataset in this redshift range. We have developed a ded- 
icated pipeline, PHASE, for extracting clean SN spectra free 
from host contamination. Our approach takes advantage of the 
rolling search mode of the SNLS by using deep stacked refer- 
ence images to estimate the host spatial profile, used during the 
SN extraction. We have also developed an identification tech- 
nique based on the simultaneous fit of light curve and spectral 
data with SALT2. We have obtained two sets of certain and prob- 
able SNe la (the SN la and SN la* categories), whose statistical 
properties have been studied in detail. We find that: 

- The statistical properties of our SN la and SN la* samples 
are similar. The SALT2 colour c and xi parameters are con- 
sistent for both subsamples, as is the "tilt" parameter, yi. 
These are indications that our samples are not strongly con- 
taminated by non SNe la, validating the inclusion of SN la* 
in the Hubble diagram. 

- The average redshift and phase for SN la* are higher than for 
SN la. In particular, the average phase of SN la* is 5.8 days, 
a phase at which the spectral S/N has decreased with respect 
to maximum and SNe la spectra closely resemble SNe Ic, 

- The uncorrected peak SN la* magnitudes appear fainter than 
the SN la, consistent with the SN la "brighter-bluer" rela- 
tionship, 

- On average, spectra of SNe identified as SN la* have a 
higher host contamination than SN la spectra, 

- We find two peculiar SNe whose spectra and properties re- 
semble those of the over-luminous SN 1999aa supernova. 

In an attempt to make a first characterisation of the physical 
properties of our spectral data, we have built average spectra for 
pre-maximum, maximum-light, and post-maximum phases, us- 
ing only SN la spectra, in two redshift bins (z < 0.5 and z > 0.5) . 
A consistent comparison of these spectra to iHsiao et aTl (l2007h 



template spectra shows a remarkable concordance in the posi- 
tions and relative flux of the spectral features. Of particular in- 
terest is the UV region of the spectra. Though more differences 
are found in this region than in the optical region, the overall 
agreement is excellent. 

We have also internally compared the average SN la spec- 
tra for z < 0.5 and z > 0.5. We find evidence for differences 
in the intermediate mass element absorptions in the average op- 
tical spectrum of SNe la between z ~ 0.35 and z ~ 0.7 (the 
average redshifts of the z < 0.5 and z > 0.5 bins). The brighter 
SNe spectra have weaker absorptions of singly ionized IMEs. 
Some discrepancy arises in the UV region of the spectra. As far 
as our spectroscopic sample is concerned, we have shown that 
these spectral differences can be essentially accounted for by the 
Malmquist bias. The use of SNe la as "calibrated candles" for 
cosmological purposes is thus fully justified. 

Fina lly, we have c ompared our VLT composite spectra with 
the .Ellis et al.l (l2008h ones, both for pre- and at maximum 
phases. We find a good agreement between the two sets of spec- 
tra, which, given their fully independent processing, gives us 
confidence in the quality of both spectral sets. 

The VLT spectra presented in this paper for which sufficient 
photometric information ex ists, once me r ged w i th the Gemini 
and Keck SNL S spe ctra jHoweU et alj 120051: iBronder et all 
|2008'; 'Ems et an i2008l) . will constitute the spectroscopic sample 
for the SNLS 3rd year cosmological analyses. 
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13 



SNname 


RA (J2000) 


DEC (J2000) 


Spectrum Date 


Exp. Time (s) 


Seeing (") 


Dsep (")* 


Air Mass 


Mag h 


03Dlar 


+02:27:14.67 


-04:19:04.74 


2003-09-28 


3x750 


0.79 


1.05 


1.21 


22.9 


03Dlbf 


+02:24:02.39 


-04:55:57.11 


2003-09-27 


3x750 


0.96 


0.79 


1.09 


23.5 


03Dlbm 


+02:24:36.07 


-04:12:44.51 


2003-09-30 


3x750 


0.55 


0.29 


1.08 


23.4 


OBDlbp 


+02:26:37.74 


-04:50:19.61 


2003-09-28 


3x750 


0.82 


1.15 


1.10 


21.9 


03D1CO 


+02:26:16.26 


-04:56:06.18 


2003-11-04 


5x750 


0.70 


* 


1.22 


23.2 


03Dldt 


+02:26:31.19 


-04:03:08.43 


2003-11-30 


3x750 


0.69 


* 


1.17 


22.4 


03Dlfc 


+02:25:43.60 


-04:08:38.61 


2003-12-25 


5x750 


0.74 


0.44 


1.32 


21.4 


03Dlfl 


+02:25:58.30 


-04:07:44.03 


2003-12-22 


5x750 


0.71 


0.38 


1.24 


22.6 


03Dlgt 


+02:24:56.01 


-04:07:36.81 


2004-01-21 


5x750 


0.77 


0.48 


1.36 


23.3 


03D4ag 


+22:14:45.83 


-17:44:22.76 


2003-06-30 


2x750 


0.56 


1.09 


1.03 


21.0 


03D4at 


+22:14:24.04 


-17:46:36.22 


2003-07-06 


3x750 


0.75 


0.15 


1.09 


22.8 


03D4au 


+22:16:09.93 


-18:04:39.32 


2003-07-06 


3x750 


0.69 


* 


1.02 


23.0 


03D4CX 


+22:14:33.75 


-17:35:15.04 


2003-09-03 


3x750 


0.54 


0.28 


1.03 


23.5 


03D4cx" 


+22:14:33.75 


-17:35:15.04 


2003-09-06 


3x750 


0.74 


* 


2.15 


23.5 


03D4cy'' 


+22:13:40.42 


-17:40:53.91 


2003-09-27 


3x750 


0.88 


1.15 


1.01 


23.6 


03D4di 


+22:14:10.24 


-17:30:23.82 


2003-09-03 


3x750 


0.59 


0.52 


1.11 


23.4 


03D4di'^ 


+22:14:10.24 


-17:30:23.82 


2003-09-07 


3x750 


0.55 


0.46 


1.04 


23.4 


03D4dy 


+22:14:50.53 


-17:57:23.32 


2003-09-30 


3x750 


0.76 


* 


1.02 


22.4 


03D4gf 


+22:14:22.93 


-17:44:02.52 


2003-11-25 


5x750 


0.70 


1.31 


1.29 


22.5 


03D4gg 


+22:16:40.19 


-18:09:52.05 


2003-11-27 


5x750 


0.51 


1.10 


1.40 


22.5 


04Dlag 


+02:24:41.10 


-04:17:19.50 


2004-01-19 


3x750 


0.83 


0.68 


1.26 


22.3 


04Dlaj 


+02:25:53.96 


-04:59:40.60 


2004-01-17 


5x750 


0.98 


* 


1.35 


22.9 


04Dlaj 


+02:25:53.96 


-04:59:40.60 


2004-01-20 


5x750 


0.79 


* 


1.34 


22.9 


04Dlak 


+02:27:33.36 


-04:19:38.78 


2004-01-24 


5x750 


0.77 


0.60 


1.38 


23.1 


04Dldc 


+02:26:18.46 


-04:18:43.09 


2004-08-12 


2x750 


0.75 


1.18 


1.23 


20.7 


04Dlff 


+02:25:38.62 


-04:54:09.47 


2004-09-11 


5x750 


0.83 


0.05 


1.07 


23.2 


04Dlhx^ 


+02:24:42.50 


-04:47:25.17 


2004-09-20 


5x750 


0.77 


0.52 


1.08 


22.8 


04Dliv 


+02:24:48.05 


-04:09:09.21 


2004-09-22 


5x750 


0.88 


0.69 


1.07 


23.5 


04Dljd 


+02:27:50.84 


-04:56:37.64 


2004-09-22 


5x750 


0.86 


* 


1.17 


23.4 


04Dlkj 


+02:27:52.66 


-04:10:48.91 


2004-10-12 


3x750 


0.54 


0.08 


1.09 


22.5 


04Dlks 


+02:24:09.52 


-04:58:43.40 


2004-10-11 


3x750 


0.75 


* 


1.16 


23.0 


04D1OW 


+02:26:42.70 


-04:18:22.28 


2004-11-14 


4x900 


0.65 


* 


1.08 


23.2 


04Dlpc 


+02:26:24.88 


-04:23:57.84 


2004-11-15 


4x900 


0.54 


0.47 


1.08 


23.5 


04Dlpd 


+02:27:39.75 


-04:37:52.22 


2004-11-13 


4x900 


0.66 


0.22 


1.07 


23.7 


04Dlpg 


+02:27:04.16 


-04:10:31.03 


2004-11-15 


3x750 


0.56 


* 


1.19 


22.8 


04Dlpp 


+02:25:12.48 


-04:42:05.70 


2004-11-14 


4x900 


0.56 


0.20 


1.11 


23.0 


04Dlqd 


+02:26:33.08 


-04:06:25.84 


2004-11-19 


4x900 


0.87 


* 


1.37 


23.1 


04Dlrh 


+02:27:47.15 


-04:15:13.56 


2004-12-11 


3x750 


0.77 


0.39 


1.66 


21.9 


04Dlrx 


+02:25:11.03 


-04:37:06.06 


2004-12-14 


4x900 


0.74 


1.60 


1.07 


23.7 


04Dlsa 


+02:27:56.16 


-04:10:33.85 


2004-12-17 


3x750 


0.97 


* 


1.12 


22.2 


04Dlsi 


+02:24:48.80 


-04:34:08.22 


2004-12-16 


4x900 


1.17 


* 


1.07 


24.0 


04D2ac 


+ 10:00:18.88 


+02:41:21.34 


2004-01-31 


3x750 


0.86 


0.29 


1.13 


21.6 


04D2al 


+ 10:01:52.48 


+02:09:50.78 


2004-01-26 


6x750 


0.83 


* 


1.26 


23.2 


04D2an 


+ 10:00:52.31 


+02:02:28.42 


2004-01-22 


5x750 


0.65 


* 


1.37 


22.7 


04D2an 


+10:00:52.31 


+02:02:28.42 


2004-01-28 


5x750 


0.70 


* 


1.14 


22.7 


04D2bt 


+09:59:32.69 


+02:14:52.99 


2004-03-20 


5x750 


0.62 


0.35 


1.19 


21.0 


04D2ca 


+ 10:01:20.49 


+02:20:21.67 


2004-03-23 


5x750 


0.72 


0.98 


1.32 


23.9 


04D2ca 


+ 10:01:20.49 


+02:20:21.67 


2004-03-21 


2x750 


0.80 


0.98 


1.48 


23.9 


04D2CC 


+09:58:44.69 


+02:16:14.53 


2004-03-22 


5x750 


0.82 


* 


1.40 


23.3 


04D2cc 


+09:58:44.69 


+02:16:14.53 


2004-03-24 


5x750 


0.74 


* 


2.15 


23.3 
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DFC f 72000^ 


vj iJt'm mil uav^ 


P xn XitTif* ( ^ 

l-jAlJ, lllllt' ) 




Dsen fl* 


All" IV^a^m 




04D2cf 


+10:01:56.09 


+01:52:46.12 


2004-03-23 


5x750 


0.79 


2.01 


1.13 


21.8 


04D2CW 


+ 10:01:22.78 


+02:11:54.81 


2004-03-21 


5x750 


0.70 


* 


1.13 


23.7 


04D2fp 


+09:59:28.12 


+02:19:15.65 


2004-04-15 


3x750 


0.66 


1.26 


1.23 


21.9 


04D2fs 


+10:00:22.09 


+01:45:55.43 


2004-04-15 


3x750 


0.60 


* 


1.13 


21.7 


04D2gc 


+10:01:39.27 


+01:52:58.99 


2004-04-15 


5x750 


0.62 


0.25 


1.13 


22.6 


04D2gp 


+09:59:20.37 


+02:30:31.96 


2004-04-20 


5x750 


0.77 


1.16 


1.20 


23.3 


04D2iu 


+ 10:01:13.18 


+02:24:54.18 


2004-05-12 


4x750 


0.92 


* 


1.48 


23.3 


04D2iu 


+ 10:01:13.18 


+02:24:54.18 


2004-05-12 


4x750 


0.78 


* 


1.14 


23.3 


04D2ja 


+09:58:48.50 


+01:46:18.25 


2004-05-13 


3x750 


0.55 


0.31 


1.15 


23.1 


04D2ja 


+09:58:48.50 


+01:46:18.25 


2004-05-14 


1x750 


0.72 


0.31 


1.12 


23.1 


04D2mc 


+ 10:00:35.56 


+01:47:03.05 


2005-01-06 


3x750 


0.84 


0.50 


1.12 


22.0 


04D4an 


+22:15:57.10 


-17:41:43.81 


2004-07-12 


5x750 


0.87 


0.29 


1.01 


23.3 


04D4bk 


+22:15:07.70 


-18:03:36.99 


2004-07-12 


5x750 


1.01 


* 


1.15 


23.3 


04D4bq 


+22:14:49.42 


-17:49:39.34 


2004-07-12 


3x750 


1.36 


1.41 


1.05 


22.5 


04D4bq 


+22:14:49.42 


-17:49:39.34 


2004-07-16 


5x750 


0.83 


1.41 


1.08 


22.5 


04D4dw 


+22:16:44.67 


-17:50:02.66 


2004-07-20 


5x750 


0.65 


* 


1.06 


23.6 


04D4fx 


+22:16:38.14 


-18:03:58.89 


2004-08-12 


5x750 


0.70 


* 


1.08 


22.7 


04D4gz 


+22:16:59.01 


-17:37:18.85 


2004-09-12 


6x750 


1.06 


0.27 


1.09 


22.7 


04D4gz 


+22:16:59.01 


-17:37:18.85 


2004-08-21 


5x750 


1.00 


0.27 


1.32 


22.7 


04D4hf 


+22:16:57.90 


-17:41:13.29 


2004-08-21 


5x750 


0.93 


* 


1.08 


23.8 


04D4ht 


+22:14:33.28 


-17:21:31.21 


2004-09-20 


5x750 


0.85 


0.30 


1.46 


21.5 


04D4ib 


+22:16:41.71 


-18:06:18.32 


2004-09-20 


5x750 


0.84 


* 


1.24 


22.7 


04D4id 


+22:16:21.40 


-17:13:44.50 


2004-09-20 


5x750 


0.85 


1.17 


1.11 


23.0 


04D4jr 


+22:14:14.32 


-17:21:00.78 


2004-10-12 


3x750 


0.64 


* 


1.01 


22.0 


04D4ju 


+22:17:02.73 


-17:19:58.13 


2004-10-12 


3x750 


0.64 


0.34 


1.04 


23.0 


04D4jw 


+22:17:18.89 


-17:39:55.78 


2004-10-12 


5x900 


0.91 


* 


1.04 


23.7 


OSDlcb 


+02:26:57.06 


-04:07:03.24 


2005-09-12 


4x900 


0.87 


1.31 


1.07 


22.8 


OSDlck 


+02:24:24.85 


-04:42:46.08 


2005-09-11 


4x900 


1.04 


0.32 


1.13 


23.3 


OSDlcl 


+02:25:05.66 


-04:13:14.95 


2005-09-27 


6x900 


0.83 


* 


1.15 


23.2 


OSDldn 


+02:24:26.65 


-04:59:29.61 


2005-09-29 


4x900 


0.81 


1.11 


1.11 


22.5 


OSDlhk 


+02:24:39.18 


-04:38:03.31 


2005-12-09 


3x750 


0.68 


* 


1.07 


20.9 


OSDlhn 


+02:24:36.24 


-04:10:54.78 


2005-12-09 


3x750 


0.85 


0.17 


1.09 


20.6 


05Dliz 


+02:26:12.77 


-04:46:24.13 


2005-12-05 


4x900 


0.71 


* 


1.08 


23.6 


OSDlke 


+02:26:46.29 


-04:51:44.99 


2005-12-30 


3x750 


0.65 


0.42 


1.09 


22.7 


05D2ac 


+09:58:59.24 


+02:29:22.21 


2005-01-08 


3x750 


0.56 


0.47 


1.16 


22.0 


05D2ay 


+ 10:01:08.32 


+02:07:49.16 


2005-01-16 


4x900 


0.70 


2.53 


1.14 


23.7 


05D2bt 


+10:01:40.24 


+02:33:58.18 


2005-02-03 


3x500 


0.64 


0.22 


1.13 


22.8 


05D2bv 


+ 10:02:17.00 


+02:14:26.16 


2005-02-03 


3x750 


0.71 


0.39 


1.14 


22.2 


05D2bw 


+09:58:30.56 


+02:32:11.04 


2005-02-03 


4x900 


0.69 


* 


1.27 


23.5 


05D2by 


+ 10:00:28.03 


+02:31:17.54 


2005-02-14 


4x900 


0.95 


* 


1.66 


23.5 


05D2cb 


+09:59:24.58 


+02:19:41.45 


2005-02-14 


2x700 


1.02 


* 


2.56 


22.6 


05D2ci 


+10:00:02.04 


+02:12:57.07 


2005-03-07 


2x900 


0.40 


0.52 


1.31 


23.0 


05D2ci 


+ 10:00:02.04 


+02:12:57.07 


2005-03-10 


2x900 


0.72 


0.52 


1.12 


23.0 


05D2ci 


+ 10:00:02.04 


+02:12:57.07 


2005-03-06 


2x900 


1.01 


0.52 


1.24 


23.0 


05D2ct 


+ 10:01:42.18 


+02:07:25.66 


2005-03-10 


8x900 


0.74 


* 


1.23 


23.6 


05D2dt 


+ 10:01:23.93 


+01:51:27.87 


2005-03-18 


4x900 


0.82 


* 


1.32 


22.9 


05D2dw 


+09:58:32.07 


+02:01:56.07 


2005-03-18 


4x900 


0.66 


1.18 


1.13 


21.8 


05D2dy 


+ 10:00:58.10 


+02:10:59.16 


2005-03-16 


4x900 


0.75 


1.29 


1.19 


22.2 


05D2eb 


+ 10:00:14.56 


+02:24:26.78 


2005-03-16 


4x900 


0.77 


0.78 


1.13 


22.4 


05D2ec 


+09:59:26.18 


+02:00:49.06 


2005-03-16 


4x900 


0.78 


0.40 


1.23 


22.8 
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Table 1 (cont'd) 



SN name 


RA (J2000) 


DEC (J2000) 


Spectrum Date 


Exp. Time (s) 


Seeins f") 


Dsep (")* 


Air Mass 


Mas im 


05D2ei 


+ 10:01:39.12 


+01:49:11.83 


2005-04-13 


3x750 


0.70 


0.63 


1.95 


21.3 


05D2fq 


+09:59:08.50 


+02:36:07.20 


2005-04-07 


4x900 


0.77 


1.26 


1.14 


23.0 


05D2he 


+ 10:01:26.51 


+02:19:03.82 


2005-04-13 


4x900 


0.83 


0.34 


1.16 


23.0 


05D2ie 


+ 10:01:02.90 


+02:39:29.11 


2005-04-09 


4x900 


0.79 


0.28 


1.33 


22.2 


05D2nn 


+09:58:49.75 


+02:42:36.95 


2005-12-07 


4x900 


0.81 


* 


1.30 


23.5 


05D4af 


+22:16:33.14 


-18:00:17.72 


2005-07-03 


4x900 


0.60 


1.49 


1.02 


22.4 


05D4ag 


+22:13:47.43 


-18:09:54.95 


2005-07-11 


4x900 


0.58 


* 


1.16 


23.1 


05D4ay 


+22:14:33.19 


-17:46:03.62 


2005-07-03 


3x900 


0.69 


* 


1.12 


23.3 


05D4ay 


+22:14:33.19 


-17:46:03.62 


2005-07-05 


1x900 


0.74 


* 


1.18 


23.3 


05D4be 


+22:16:53.39 


-17:14:09.96 


2005-07-07 


3x750 


0.56 


0.18 


1.03 


22.1 


05D4bi 


+22:15:56.55 


-17:59:09.85 


2005-07-07 


4x900 


0.64 


* 


1.11 


22.9 


05D4bj 


+22:14:59.66 


-18:08:00.26 


2005-07-10 


4x900 


0.65 


* 


1.02 


23.0 


05D4cn 


+22:13:31.45 


-17:17:19.84 


2005-08-10 


4x900 


0.74 


1.00 


1.10 


23.0 


05D4cq 


+22:14:09.65 


-18:13:36.17 


2005-08-04 


3x750 


0.57 


1.29 


1.01 


22.8 


05D4CS 


+22:17:17.82 


-17:52:48.54 


2005-08-04 


3x750 


0.90 


* 


1.02 


22.8 


05D4CW 


+22:14:50.08 


-17:44:19.56 


2005-08-08 


3x750 


0.64 


* 


1.03 


21.6 


05D4dw 


+22:16:55.88 


-18:03:05.26 


2005-09-10 


6x900 


0.83 


0.28 


1.05 


23.4 


05D4ef 


+22:13:59.31 


-18:12:44.44 


2005-09-24 


6x900 


1.53 


* 


1.47 


23.0 


05D4ej 


+22:15:52.52 


-18:11:44.59 


2005-09-27 


4x900 


0.65 


0.03 


1.05 


22.9 


05D4ek 


+22:16:27.52 


-17:44:10.74 


2005-09-26 


4x900 


0.77 


1.19 


1.06 


22.6 


05D4ev 


+22:13:30.42 


-17:42:02.83 


2005-09-12 


6x900 


0.49 


2.87 


1.02 


23.4 


05D4fe 


+22:15:40.27 


-17:31:42.38 


2005-09-27 


4x900 


0.60 


* 


1.01 


23.9 


05D4ff 


+22:16:20.16 


-18:02:33.20 


2005-09-29 


4x900 


0.79 


* 


1.12 


22.0 


05D4fg 


+22:16:41.35 


-17:35:43.66 


2005-09-29 


4x900 


0.74 


0.56 


1.09 


23.0 


06Dlab 


+02:24:59.03 


-04:40:51.20 


2006-01-06 


3x750 


1.05 


0.53 


1.15 


21.5 


06D2bk 


+09:58:42.87 


+02:10:19.14 


2006-02-09 


3x750 


0.72 


0.41 


1.18 


22.6 


06D2ca 


+09:59:09.56 


+02:18:09.06 


2006-02-27 


A . .r\f\f\ 

4x900 


0.91 


0.33 


1.27 


22.5 


06D2ca 


+09:59:09.56 


+02:18:09.06 


2006-02-28 


1x900 


0.90 


0.33 


1.32 


22.5 


06D2cb-'^ 


+ 10:01:23.49 


+02:27:22.36 


2006-03-02 


6x900 


0.88 




1.13 


24.1 


06D2cc 


+09:59:26.89 


+02:27:04.88 


2006-02-28 


4x900 


0.80 


1.42 


1.18 


22.6 


06D2cd"'^ 


+ 10:00:57.96 


+02:32:02.54 


2006-03-04 


6x900 


0.60 


* 


1.24 


24.0 


06D2ce'^'^ 


+ 10:01:44.22 


+02:42:47.93 


2006-03-03 


5x900 


0.76 


* 


1.26 


23.3 


06D2ck 


+ 10:01:28.02 


+01:51:49.07 


2006-03-08 


4x900 


1.07 


0.34 


1.30 


22.7 


06D2ga 


+09:59:08.11 


+01:56:18.58 


2006-05-25 


4x900 


0.66 


3.60 


1.27 


23.4 


06D4ce 


+22:17:17.25 


-17:14:09.37 


2006-07-19 


4x900 


0.86 


* 


1.23 


23.1 


06D4ce 


+22:17:17.25 


-17:14:09.37 


2006-07-20 


4x900 


0.79 


* 


1.04 


23.2 


06D4cl 


+22:16:57.05 


-17:27:02.28 


2006-07-21 


8x900 


0.70 


4.34 


1.03 


23.4 


06D4CO 


+22:15:26.52 


-17:52:09.07 


2006-07-30 


3x750 


1.08 


* 


1.01 


21.8 


06D4cq 


+22:16:55.52 


-17:42:43.13 


2006-07-30 


4x750 


1.07 


0.84 


1.10 


21.8 



" Except for 8 SN la taken in MOS mode. 

* Host-supernova distance (from PHASE model) in arcsec. An asterisk either denotes cases for which a separate extraction of the 
supernova and the host component was not possible with PHASE or the absence of a detectable host. 
Observed with Grism 3001. 

SNLS-03D4cy is identified as a SN la on the basis of a Gemini spectrum published in lHowell et al.l (l2005l) . The VLT spectrum is 
shown here for reference. 

There are two underlying galaxies in the line-of-sight of SNLS-04Dlhx. 
f FORS2 spectrum 



16 Balland et al.: The ESOA'^LT 3rd year Type la supernova data set from the SNLS 

Table 2. Results of identifications of the SNLS 3rd year ESO-VLT SN la and SN la* spectra. 



SNname 


Type 


z 


z source 


0) 


Host model" 


Host fraction 


<S/N: 


03Dlar 


la* 


0.408±0.001 


H 


5.3 


Ed) 


0.16 


2.3 


03Dlbf 


la* 


703+0 001 


H 


-2.7 




0.16 


1.4 


OSDlbtn 


la* 


iTiiO 001 


H 


-5.1 




0.07 


2.1 


OSDIbn 


la* 


347+0 001 


H 


-7.5 




0.86 


2.2 


03Dlco 


la 


0.679±0.001 


H 


-4.1 




0.10 


2.0 


03Dldt 


la 


612+0 001 


H 


5.1 




0.65 


3.3 


03Dlfc 


la 


332+0 001 


H 


-4.4 




0.33 


21.7 


03Dlfl 


la 


687±0 001 


H 


0.5 






3.0 


03Dlgt 


la 


560+0 001 


H 


7.0 




0.50 


2.4 


03D4ag 


To 

^"pec 


0.285+0.001 


H 


-8.6 


E(2) 


0.27 


29.5 


03D4at 


la 


634+0 001 


H 


5.5 


Sb-Sc 


0.42 


3.7 


03D4au 


la* 


468±0 001 


H 


6.5 


Sb-Sc 


0.43 


5.2 


03D4cx 


la 


0.949+0.014 


s 


1.2 


NoGalaxy 




0.6 


03D4cx 


la 


0.949+0.014 


s 


2.7 


S0(4) 


0.95 


1.5 


03D4cy 


la 


0.927+0.001 


H 


4.6 


Sb-Sc 


0.16 


0.5 


03D4di 


la* 


0.899+0.001 


H 


-8.6 


Sb-Sc 


0.11 


1.8 


03D4di 


la* 


0.899+0.001 


H 


-6.5 


NoGalaxy 




1.8 


03D4dy 


la 


61+0 01 


s 


4.8 


IVodalaxv 




5.5 


03D4gf 


la* 


0.580+0.001 


H 


20.4 




0.46 


1.5 




la* 


0.592+0.001 


H 


16.4 


Sc(lO) 


0.73 


1.7 


04Dlag 


la 


0.557+0.001 


H 


4.3 


EH 2) 


0.08 


5.4 


04Dlaj 


la* 


0.721+0.011 


s 


11.6 


NoGalaxy 




1.0 


04Dlaj 


la* 


0.721+0.011 


s 


13.4 


NoGalaxy 




1.2 


04Dlak 


la* 


0.526+0.001 


H 


11.8 


NoGalaxy 




0.6 


04Dldc 


la 


0.211+0.001 


H 


-0.4 


NoGalaxy 




33.4 


04DlfF 


la 


0.86+0.01 


s 


4.6 


Sc(2) 


0.24 


2.0 


04Dlhx 


la 


560+0 001 


H 


5.9 


NoGalaxv 




1.9 


04Dliv 


la 


0.998+0.001 


H 


3.0 


NoGalaxy 




1.3 


04Dljd 


la* 


0.778+0.001 


H 


6.8 


Ef2) 


0.58 


1.9 


04Dlkj 


la 


0.585+0.001 


H 


-3.7 


NoGa.la.xy 




6.7 


04Dlks 


la 


0.798+0.012 


s 


-1.0 


NoGalaxy 




2.8 


04D1ow 

\j r \^ J. V/ vv 


la 


915+0 001 


H 


6.4 


NnCralflxv 




1.6 


04Dlnc 


la 


770+0 001 

\j . / / \j -J—\J .vjyji. 


H 


0.1 


EH) 


0.22 


3.0 


04DlDd 


la 


95+0 01 


s 


2.5 


Ed) 


0.31 


2.1 




la 


515+0 001 


H 


-1.3 


Sa-Sb 


0.61 


8.9 


04Dlpp 


la 


0.735+0.001 


H 


2.2 


NoGalaxy 




1.2 


04Dlqd 


la* 


Q 767+0 001 


H 


-0.2 


ECl) 


0.71 


4.8 


04Dlrh 


la 


Q 436+0 001 

\J .~^\J -!—\J .XJXJ 1. 


H 


0.0 


NoGalaxv 

L y V 1 V.J cixciA y 




7.8 


04Dlrx 


la* 


985+0 001 


H 


0.9 


NoGalaxv 




1.4 


04Dlsa 


la 


0.585+0.001 


H 


-2.6 


EM) 


0.72 


6.4 


04Dlsi 


la 


0.702+0.001 


H 


-1.7 


Sa('9) 


0.79 


7.6 


04D2ac 


la 


348+0 001 

\J . •^~\J -!—\J .\J\J 1. 


H 


1.3 


Sa-Sb 


0.20 


7.6 


04D2al 


la 


0.836+0.001 


H 


-2.5 


E(2) 


0.29 


2.5 


0407 an 


Ta 


62+0 01 




-3 4 


^Tor^fiiiiYv 




7 7 


04D2an 


la 


0.62+0.01 


S 


0.3 


Sc(9) 


0.06 


7.4 


04D2bt 


la 


0.220+0.001 


H 


6.6 


E(6) 


0.01 


9.2 


04D2ca 


la* 


0.835+0.001 


H 


12.0 


NoGalaxy 




0.5 


04D2ca 


la* 


0.835+0.001 


H 


13.1 


NoGalaxy 




0.9 


04D2CC 


la* 


0.838+0.001 


H 


6.1 


S0(3) 


0.66 


1.8 


04D2CC 


la* 


0.838+0.001 


H 


7.2 


80(2) 


0.67 


1.5 


04D2cf 


la 


0.369+0.001 


H 


8.5 


NoGalaxy 




12.1 


04D2CW 


la* 


0.568+0.001 


H 


19.3 


Sb-Sc 


0.58 


1.4 
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Table 2 (cont'd) 



SN name Type z z source 



04D2fp 


la 


0.415±0.001 


H 


U4Uzrs 


la 




H 


04D2gc 


la 


0.521±0.001 


T T 

H 


04D2gp 


la 


0.732±0.001 


H 


04D2iu 


la* 


0.70±0.01 


S 


04D2iu 


la* 


0.70±0.01 


S 


04D2ja 


la* 


0.740+0.001 


H 


04D2ja 


la* 


0.740±0.001 


H 


04D2mc 


la 


f\ ^ A n ^ /A /A/A -1 

0.348±0.001 


H 


04D4an 


la 


/A -i -"^ /A /A/A-1 

0.613+0.001 


H 


04D4bk 


la* 


/A O O /A /A 1 

0.88+0.01 


S 


04D4Dq 


la 


0.55+0.01 


s 


04D4bq 


la 


A C C ■ A A 1 

0.55±0.01 


s 


04D4dw 


la* 


-1 /AO -1 . /A /A/A-1 

1.031+0.001 


H 


04D4fx 


la 


/A y'/A/A /A /A/A-1 

0.629+0.001 


H 


04D4gz 


la* 


r\ '~i C . f\ /A/A 1 

0.375+0.001 


H 


04D4gz 


la* 


/A ^ ^ C /A /A /A -1 

0.375+0.001 


H 


04D4hf 


la* 


/A y . f\ /A -1 A 

0.936±0.014 


S 


f\ A -l-v At 1 

04D4ht 


la 


/A -1 ^7 . /A /A/A 1 

0.217+0.001 


H 


04D4ib 


la 


/A ^/A/A /A /A/A1 

0.699+0.001 


H 


04D4id 


la 


0.769+0.001 


H 


04D4jr 


la 


0.47+0.01 


S 


04D4ju 


la 


/A >l . /A ^/A -1 

0.472+0.001 


H 


U4U4JW 


la* 


r\ r\£ii 1 r\ r\Ai 


T T 

H 


05Dlcb 


la 


t\ ^ , r\ r\r\-\ 

0.632+0.001 


H 


OSDlck 


la 


^ 1 ^ . /A /A /A 1 

0.617+0.001 


H 


OSDlcl 


la* 


/A O -^ . /A /A -1 

0.83+0.01 


S 


05Dldn 


la 


A C£l£l 1 A AA1 

0.566+0.001 


TT 


OSDlhk 


lapec 


/A /A^^O . /A ^/A 1 

0.263+0.001 


H 


05Dlhn 


la 


1 -1 /A . /A /A/A 1 

0.149+0.001 


H 


OSDliz 


la* 


0.86+0.01 


S 


OSDlke 


la 


/A ✓'/A . /A /A -1 

0.69+0.01 


s 


05D2ac 


la 


0.479+0.001 


H 


05D2ay 


la* 


r\ /A'A . /A /A 1 

0.92+0.01 


S 


05D2bt 


la 


/A O /A/A-1 

0.68+0.01 


s 


UjUzdv 


la 


A AHA 1 A AA1 


1 1 

H 


05D2bw 


la 


/A JA/^ . /A /A -1 

0.92+0.01 


s 


05D2by 


la* 


/A O /A -i /A /A /A -1 

0.891+0.001 


H 


05D2cb 


la* 


/A -1 . /A /A/A 1 

0.427+0.001 


H 


U5U2C1 


la 


A O A 1 A AA 1 

0.630+0.001 


T T 

H 


05D2ci 


la 


0.630+0.001 


H 


05D2ci 


la 


0.630+0.001 


H 


05D2ct 


la* 


0.734+0.001 


H 


05D2dt 


la 


0.574+0.001 


H 


05D2dw 


la 


0.417+0.001 


H 


05D2dy 


la 


0.51+0.01 


S 


05D2eb 


la 


0.534+0.001 


H 


05D2ec 


la 


0.640+0.001 


H 


05D2ei 


la* 


0.366+0.001 


H 


05D2fq 


la 


0.733+0.001 


H 


05D2he 


la 


0.608+0.001 


H 



VP 


Host model" 


Host fraction 


^ C / AT" ■ 
< o /iv . 


1.8 


NoGalaxy 


- 


11.5 


1.7 


NoGalaxy 


- 


17.8 


-4.9 


E(l) 


0.19 


5.2 


2.7 


NoGalaxy 


- 


2.5 


9.8 


NoGalaxy 


- 


0.7 


10.4 


E(l) 


0.09 


1.2 


9.2 


E(3) 


0.44 


0.8 


9.8 


E(3) 


0.19 


0.5 


6.5 


Sa-Sb 


0.30 


10.3 


8.1 


NoGalaxy 


- 


0.8 


3.1 


E(l) 


0.20 


1.6 


2.6 


NoGalaxy 


- 


2.6 


5.1 


NoGalaxy 


- 


6.5 


2.1 


Sa-Sb 


0.68 


1.5 


-8.4 


NoGalaxy 


- 


5.8 


-5.8 


Sa-Sb 


0.33 


2.2 


10.2 


NoGalaxy 


- 


1.8 


-0.5 


NoGalaxy 


- 


0.8 


6.6 


Sa-Sb 


0.62 


3.9 


0.7 


SO-Sa 


0.71 


4.2 


3.0 


NoGalaxy 


- 


2.3 


-5.9 


NoGalaxy 


- 


16.5 


-2.2 


S0(1) 


0.44 


3.5 


2.2 


SO-Sa 


0.69 


2.3 


4.3 


Sb-Sc 


0.28 


4.5 


-2.7 


Sb-Sc 


0.39 


0.7 


8.7 


Sd(4) 


0.36 


1.8 


-4.7 


NoGalaxy 


- 


7.3 


-5.0 


Sc(5) 


0.39 


21.4 


-1.0 


NoGalaxy 


- 


13.0 


7.8 


Sd(l) 


0.24 


2.1 


2.2 


NoGalaxy 


- 


4.9 


2.0 


NoGalaxy 


- 


10.9 


5.8 


S0(1) 


0.07 


1.1 


1.4 


NoGalaxy 


- 


2.1 


-0.1 


Sb-Sc 


0.26 


8.3 


1.9 


Sb-Sc 


0.25 


2.1 


0.9 


Sc(ll) 


0.83 


1.3 


-5.8 


Sc(ll) 


0.57 


2.9 


3.6 


NoGalaxy 


- 


0.8 


4.9 


NoGalaxy 


- 


0.7 


6.7 


NoGalaxy 


- 


1.4 


7.8 


Sd(ll) 


0.63 


4.6 


-1.7 


E(3) 


0.43 


6.0 


-5.3 


S0(1) 


0.12 


10.3 


1.1 


Sb-Sc 


0.02 


4.8 


-4.7 


S0(1) 


0.04 


6.2 


2.7 


NoGalaxy 




2.4 


16.9 


S0(4) 


0.55 


2.5 


2.1 


NoGalaxy 




1.9 


3.0 


E(l) 


0.26 


3.1 
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Table 2 (cont'd) 



SN name 


Type 


z 


z source 


O 


Host model" 


Host fraction 


< SIN 


05D2ie 


la 


0.348+0.001 


H 


-8.9 


Sc(4) 


0.45 


4.4 


05D2nn 


la 


0.87+0.01 


S 


-1.3 


Sa(8) 


0.26 


1.9 


05D4af 


la 


0.499+0.001 


H 


11.3 


NoGalaxy 


- 


8.8 


05D4ag 


la* 


0.64+0.01 


S 


14.0 


Sb-Sc 


0.10 


2.0 


05D4ay 


la* 


0.408+0.001 


H 


8.2 


Sa-Sb 


0.80 


2.55 


05D4ay 


la* 


0.408+0.001 


H 


9.6 


Sa-Sb 


0.80 


1.8 


05D4be 


la 


0.537+0.001 


H 


3.9 


E(l) 


0.31 


5.5 


05D4bi 


la 


0.775+0.001 


H 


-1.5 


Sc(l) 


0.40 


3.7 


05D4bj 


la 


0.701+0.001 


H 


-2.0 


Sc(9) 


0.68 


2.8 


05D4cn 


la 


0.763+0.001 


H 


4.7 


S0(1) 


0.27 


3.0 


05D4cq 


la 


0.701+0.001 


H 


-0.1 


NoGalaxy 


- 


4.6 


05D4CS 


la 


0.79+0.01 


S 


-1.7 


Sc(12) 


0.16 


4.3 


05D4CW 


la 


0.375+0.001 


H 


7.0 


S0(6) 


0.57 


16.4 


05D4dw 


la 


0.855+0.001 


H 


4.7 


Sc(12) 


0.19 


0.9 


05D4ef 


la* 


0.605+0.001 


H 


3.5 


E(3) 


0.78 


1.6 


05D4ej 


la 


0.585+0.001 


H 


7.6 


NoGalaxy 


- 


3.5 


05D4ek 


la 


0.536+0.001 


H 


2.1 


NoGalaxy 


- 


8.2 


05D4ev 


la 


0.722+0.001 


H 


-3.6 


NoGalaxy 


- 


3.4 


05D4fe 


la* 


0.984+0.001 


H 


-2.0 


E(3) 


0.65 


1.1 


05D4ff 


la* 


0.402+0.001 


H 


5.1 


Sc(6) 


0.83 


4.9 


05D4fg 


la 


0.839+0.001 


H 


-0.3 


S0(5) 


0.32 


2.3 


06Dlab 


la 


0.182+0.001 


H 


-4.4 


E(7) 


0.11 


13.9 


06D2bk 


la 


0.499+0.001 


H 


0.9 


Sb-Sc 


0.15 


3.8 


06D2ca 


la 


0.531+0.001 


H 


0.0 


Sd(5) 


0.27 


3.2 


06D2ca 


la 


0.531+0.001 


H 


0.7 


Sd(3) 


0.18 


2.5 


06D2cb 


la* 


1.00+0.01 


S 


3.0 


NoGalaxy 


- 


1.7 


06D2CC 


la 


0.532+0.001 


H 


3.3 


NoGalaxy 


- 


6.2 


06D2cd 


la 


0.930+0.001 


H 


4.2 


NoGalaxy 


- 


5.8 


06D2ce 


la 


0.82+0.01 


S 


0.0 


NoGalaxy 




7.3 


06D2ck 


la* 


0.552+0.001 


H 


7.7 


Sa-Sb 


0.24 


1.6 


06D2ga 


la* 


0.84+0.01 


S 


5.4 


Sb-Sc 


0.19 


2.0 


06D4ce 


la 


0.85+0.01 


S 


2.7 


NoGalaxy 




2.3 


06D4ce 


la 


0.85+0.01 


s 


3.3 


NoGalaxy 




2.7 


06D4cl 


la 


1.00+0.01 


s 


-2.5 


E(l) 


0.01 


2.5 


06D4CO 


la 


0.437+0.001 


H 


3.5 


E-SO 


0.58 


12.8 


06D4cq 


la 


0.411+0.001 


H 


-1.4 


NoGalaxy 




8.2 



° Host type of the best fitting model. The number in parenthesis is the age of the best fitting PEGA SE template (in Gyrs ). Two host 
types separated by a dash indicate that the best fitting host model is an interpolation between two iKinney et al.l (Il996h templates. 
The label 'NoGalaxy' denotes that the best fitting model is obtained with no host model. 
Average signal-to-noise ratio per 5 A pixel. 
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Table 3. Distribution of t5?pes per SNLS Deep field and in total 
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Dl 


D2 


D3 


D4 


Total number of SNe (spectra) 


SN la 


27 


29 





30 


86 (93) 


SNIa* 


11 


13 





14 


38 (46) 
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Table 4. Average properties of the SN la and SN la* subsamples 





<z> (o-^) 


<<f>> icr<p) 


<C> (CTc) 


< ^1 > (0-;c,) 


SN la 
SNIa* 


0.60+0.02 (0.20) 
0.70±0.03 (0.21) 


1.1+0.4(4.1) 
5.8±1.1 (7.1) 


0.009+0.013(0.12) 
0.05±0.03 (0.16) 


-0.02+0.09 (0.9) 
0.004±0.18(1.1) 




<ri > (o"ri) 




fgali(^fsai) 


<5/A^>'' 




0.26+0.11 (1.1) 

0.26+0.29(1.9) 


23.95+0.04(0.36)* 

24.12+0.08(0.50)'' 


0.32+0.03 (0.23) 

0.47+0.05 (0.28) 


5.8+0.6 

1.8+0.2 


1.13 (0.30)*^ 

1.03 (0.16)'' 



" Average signal-to-noise ratio per 5 A pixel. 

* Distance corrected magnitudes in the B band, using Qq = 0.27, Qa = 0.73 and Hq = 70 km/s/Mpc. 
Excluding > 2 events (1.24 (0.49) with all events). 
Excluding > 2 events (1.46 (2.30) with all events). 
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Table 5. Properties of SNLS-03D4ag and SNLS-05Dlhk VLT spectra 



SN name 


z 




c 


xi(Am*5) 


n 


03D4ag 


0.285 


-9 


-0.057 


0.808(0.89) 


-0.02 


05Dlhk 


0.263 


-5 


-0.01 


1.367(0.97) 


-1.25 



" Phase in days with respect to B-ban d maximum light. 
Aoti5 value computed from x\ using iGuv etalJ(l2007h formula. 
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Table 6. Number of SN la spectra in each phase bin (in days) used to create VLT composite spectra. 



< -4 -4< <4 > 4 Total 



z<0.5 


7 


12 


6 


25 


z>0.5 


5 


45 


18 


68 


Total 


12 


57 


24 


93 
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Fig. 1. Redshift distribution of SNLS candidates identified as SN la or SN la*. In light green, the histogram for all SNe la 
(SN la+SN la*). In dark green, the histogram for SNe la alone. The average redshifts for the SNe la and SNe la* subsamples 
are < z >sNia- 0.60 + 0.02 and < z >sN!a*- 0.70 + 0.03 for 86 and 38 SNe, respectively. 




Fig. 2. Phase distribution of SNLS candidates identified as SN la or SN la*. In light green, the histogram for all SNe la 
(SN la+SN la*). In dark green, the histogram for SNe la alone. The average phases for the SNe la and SNe la* subsamples 
are < >sn la- 1-1 ± 0.4 days and < (f> >sn la*- 5.8 + 1.1 days, respectively. 
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Fig. 3. Colour distribution of SNLS candidates identified as SN la or SN la*. In light green, the histogram for all SNe la 
(SN la+SN la*). In dark green, the histogram for SNe la alone. The average colours for the SNe la and SNe la* subsamples 
are < c >sNia= 0.009 ± 0.013 mag and < c >sNia*= 0.05 + 0.03 mag, respectively. 




Fig. 4. SALT2 xi distribution of SNLS candidates identified as SN la or SN la*. In light green, the histogram for all SNe la 
(SN la+SN la*). In dark green, the histogram for SNe la alone. The average xi for the SNe la and SNe la* subsamples are 
< xi >sNia- -0.02 + 0.09 and < xi >sNia*- 0.004 + 0.18, respectively. 
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Fig. 5. SNLS-03D4ag (top) and SNLS-05Dlhk (bottom) rest frame spectra compared to SN 1999aa templates at -8 and -5 days 
respectively (from Matheson et al. 2008). SN 1999aa spectra are not de-redshifted. None of the spectra are corrected for peculiar 
velocities, which explains the shifts observed in the positions of absorption features. In both cases, note the very blue spectrum and 
shallow Si II /I 6150. 
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Fig. 6. Illustration of the magnitude Am of the average recalibration ("tilt") needed for the photometric models to fit the spectra, as a 
function of wavelength. For the purpose of this illustration, we have selected a subset of VLT spectra with redshifts between z - 0.4 
and z - 0.6. The average redshift is then < z >- 0.512 for SN la and < z >- 0.518 for SN la*. Errors shown in the figure (in 
light blue for SN la* and in dark blue for SN la), are computed from the average value of the recalibration parameter y\, estimated 
from the subset, for both SN la and SN la* populations. It appears that, if recalibration is about the same, on average, for SN la and 
SN la* at this redshift, the yi values are more dispersed for the SN la*. The global flux correction is mild, k\5% at both ends of 
the wavelength scale. 
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Fig. 7. Distribution of S ALT2 xl values for the whole (SN la and SN la*) sample. The vertical red/grey solid line indicates the cut 
atxl = 2 used to identify potential non- or pecuhar SN la. The two highest objects are not shown. 
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Fig. 8. In grey /blue: Average spectrum built from all VLT SN la spectra of this study for z < 0.5 (left column) and z > 0.5 (right 
column), and for pre-maximum (O < -4 days, top), around maximum (-4 days < O < +4 days, middle) and past maximum 
(<£ > +4 days, bottom). To construct the average spectrum in each panel, all individual spectra have been brought into the rest 
frame and rebinned with 5 A. Spectra have been colour-corrected using the SALT2 colour law and the colour parameter c for each 
individual spectrum. Fluxes have then been normalised to the same integral around 4500 A. For each wavelength bin, an average 
flux fx and a weighted Icr error are computed from all the spectra with measured flux in this bin. Th e number and aver age phase 
thus vary at each wavelength step. In red: an average spectrum built in the same way using this time the lHsiao et al.l ( l2007h template 
series. For each wavelength bin, the same number of individual templates and same weights as in the average VLT spectrum have 
been used, with the corresponding phases rounded off to the nearest integer value. 
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Fig. 9. Comparison of z < 0.5 (in black) and z > 0.5 (in grey/blue) VLT average spectra for pre-maximum phases (top), around 
maximum phases (middle), and post-maximum phases (bottom). In each case, the 50 A binned residuals A = /'-"^(/l)-/'^"^(/l) are 
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Fig. 10. Comparison of z < 0.5 (in black) and z > 0.5 (in grey/blue) VLT average spectra around maximum light for bluer than 
average (c < 0) subsets (top panel). The 50 A b inned residuals A - /~-°^(/l) - f^^'^-^{A) are shown in the bottom panel for the VLT 
spectra (in black) and for the lHsiao et alj ( |2007| templates (in red). For the VLT residual, variances computed as the quadratic sum 
of the variances of the two average spectra (z < 0.5 and z > 0.5) and +lcr errors are shown. Compared to the middle panel of Fig.|9] 
the differences from 3500 to 4500 A between z < 0.5 and z > 0.5 have been drastically reduced. VLT spectra have been individually 
colour corrected using the SALT2 colour law, recalibrated and normalised to the same integral in the range 4450-4550A. 
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Fig. 11. Comparison of th e VLT average spe ctra obtained from a subset of 51 spectra in the range z - 0.35 - 0.7 (black), with the 
Keck average spectra froml Ellis et alJ (1200 8*) (in grey/red), for pre-maximum phases (top), and at maximum phases (bottom). VLT 
spectra have been individually colour corrected using the SALT2 colour law and recalibrated at the first order using the yi 'tilt' 
coefficient. All spectra have been normalised to unity at 3900 A. For each phase range, the residuals A{A) - A^^^(A) - A^''^'^{A) are 
shown. Variances are computed as the quadratic sum of the variances of the two average spectra (VLT and Keck). Errors are +2cr. 
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Fig. A.l. SNLS- 03Dlar270: a SN la* supernova at z=0.408. The spectrum phase is 5.3. A E(7) host model has been subtracted. 
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Fig. A.2. SNLS- 03Dlbf269: a SN la* supernova at z=0.703. The spectrum phase is -2.7. A E(3) host model has been subtracted. 
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Fig. A.3. SNLS- 03Dlbm272: a SN la* supernova at z=0.575. The spectrum phase is -5.1. A E(3) host model has been subtracted. 



Rest-frame wavelength (A) 

„ 3000 3500 4000 4500 5000 5500 6000 

2.5r-r 




Rest-frame wavelength (A) 
3000 3500 4000 4500 5000 5500 6000 



4000 4500 5000 5500 6000 6500 7000 7500 8000 8500 
Observed Wavelength (A) 




■4000 4500 5000 5500 6000 6500 7000 7500 8000 8500 
Observed Wavelength (A) 



Fig. A.4. SNLS- 03Dlbp270: a SN la* supernova at z=0.347. The spectrum phase is -7.5. A E(2) host model has been subtracted. 
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Fig. A.5. SNLS- 03Dlco307: a SN la supernova at z=0.679. The spectrum phase is -4.1. A E(8) host model has been subtracted. 
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Fig. A.6. SNLS- 03Dldt333: a SN la supernova at z=0.612. The spectrum phase is 5.1. A Sc(4) host model has been subtracted. 
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Fig. A.7. SNLS- 03Dlfc358: a SN la supernova at z=0.332. The spectrum phase is -4.4. A E(3) host model has been subtracted. 
Three re-calibration parameters have been used to fit this spectrum. 
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Fig. A.8. SNLS- 03Dlfl355: a SN la supernova at z=0.688. The spectrum phase is 0.5. Best-fit obtained for a model with no host 
galaxy component. 
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Fig. A.9. SNLS- 03Dlgt385: a SN la supernova at z=0.560. The spectrum phase is 7.0. A E(4) host model has been subtracted. 
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Fig. A.IO. SNLS- 03D4agl80: a SN Ia_pec supernova at z=0.285. The spectrum phase is -8.6. A E(2) host model has been sub- 
tracted. Note the poor fit below rest-frame 4000 A of this SN 1999aa-lik;e supernova. 
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Fig. A.ll. SNLS- 03D4atl86: a SN la supernova at z=0.634. The spectrum phase is 5.5. A Sb-Sc host model has been subtracted. 




Fig. A.12. SNLS- 03D4aul86: a SN la* supernova at z=0.468. The spectrum phase is 6.5. A Sb-Sc host model has been subtracted. 
Note the UV flux excess in the spectrum compared to the model. This is due to a host galaxy residual in this region of the spectrum. 



Balland et al.: The ESO/VLT 3rd year Type la supernova data set from the SNLS, Online Material p 9 



Rest-frame wavelength (A) 
2000 2500 3000 3500 4000 




■foOO 4500 5000 5500 6000 6500 7000 7500 8000 8500 
Observed Wavelength (A) 



Fig. A.13. SNLS- 03D4cx245: a SN la supernova at z=0.949. The spectrum phase is 1.2. Best-fit obtained for a model with no host 
galaxy component. 
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Fig. A.14, SNLS- 03D4cx248: a SN la supernova at z=0.949. Spectrum obtained with the 3001 Grism. The spectrum phase is 2.7. 
A S0(4) host model has been subtracted. 
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Fig. A.15. SNLS- 03D4cy269: a SN la supernova at z=0.927. The spectrum phase is 4.6. A Sb -Sc host model has been subtracted. 
The identification of this supernova as a SN la is based on a Gemini spectrum of this same SN (iHowell et al.ll2005h 




Fig, A.16. SNLS- 03D4di245: a SN la* supernova at z=0.899. The spectrum phase is -8.6. A Sb-Sc host model has been subtracted. 




Fig. A.17. SNLS- 03D4di249: a SN la* supernova at z=0.899. Spectrum obtained with the 3001 Grism. The spectrum phase is -6.5. 
Best-fit obtained for a model with no host galaxy component. 




Fig. A.18. SNLS- 03D4dy272: a SN la supernova at z=0.61. The spectrum phase is 4.8. Best-fit obtained for a model with no host 
galaxy component. 
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Fig. A.19. SNLS- 03D4gf329: a SN la* supernova at z=0.580. The spectrum phase is 20.4. A S0(1) host model has been subtracted. 
The UV part of the model is poorly constrained due to the lack of UV templates at late phases in the version of the SALT2 training 
set used in this work. 




Fig. A.20. SNLS- 03D4gg330: a SN la* supernova at z=0.592. The spectrum phase is 16.4. A Sc(lO) host model has been sub- 
tracted. The UV part of the model is poorly constrained due to the lack of UV templates at late phases in the version of the SALT2 
training set used in this work. 
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Fig. A.21. SNLS- 04Dlag383: a SN la supernova at z=0.557. The spectrum phase is 4.3. A E(12) host model has been subtracted. 




Fig. A.22. SNLS- 04Dlaj381: a SN la* supernova at z=0.72L The spectrum phase is n.6. Best-fit obtained for a model with no 
host galaxy component. 
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Fig. A.23. SNLS- 04Dlaj384: a SN la* supernova at z=0.721. The spectrum phase is 13.4. Best-fit obtained for a model with no 
host galaxy component. 




Fig. A.24. SNLS- 04Dlak388: a SN la* supernova at z=0.526. The spectrum phase is 1 1.8. Best-fit obtained for a model with no 
host galaxy component. 




Fig.A.25. SNLS- 04Dldc589: a SN la supernova at z=0.211. The spectrum phase is -0.4. Best-fit obtained for a model with no 
host galaxy component. 
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Fig. A.26. SNLS- 04DlfF619: a SN la supernova at z=0.86. The spectrum phase is 4.6. A Sc(2) host model has been subtracted. 




Fig. A.27. SNLS- 04Dlhx628: a SN la supernova at z=0.560. The spectrum phase is 5.9. Best-fit obtained for a model with no host 
galaxy component. Note the model is unable to reproduce the large feature at rest-frame 4500 A. 




Fig. A.28. SNLS- 04Dliv630: a SN la supernova at z=0.998. The spectrum phase is 3.0. Best-fit obtained for a model with no host 
galaxy component. 



Balland et al.: The ESOA'LT 3rd year Type la supernova data set from the SNLS, Online Material p 17 



-0.5- 



Rest-frame wavelength (A) 
3000 4000 




Rest-frame wavelength (A) 
3000 4000 



"4000 4500 5000 5500 6000 6500 7000 7500 8000 8500 
Observed Wavelength (A) 




"4000 4500 5000 5500 6000 6500 7000 7500 8000 8500 
Observed Wavelength (A) 



Fig. A.29. SNLS- 04Dljd630: a SN la* supernova at z-O.llS. The spectrum phase is 6.8. A E(2) host model has been subtracted. 




Fig. A.30. SNLS- 04Dlkj650: a SN la supernova at z=0.585. The spectrum phase is -3.7. Best-fit obtained for a model with no host 
galaxy component 
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Fig. A.31. SNLS- 04Dlks649: a SN la supernova at z=0.80. The spectrum phase is -1.0. Best-fit obtained for a model with no host 
galaxy component. 
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Fig. A.32. SNLS- 04Dlow683: a SN la supernova at z=0.915. The spectrum phase is 6.4. Best-fit obtained for a model with no 
host galaxy component. 




Fig. A.33. SNLS- 04Dlpc684: a SN la supernova at z-Q.llQ. The spectrum phase is 0.1. A E(3) host model has been subtracted. 
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Fig. A.34. SNLS- 04Dlpd682: a SN la supernova at z=0.95. The spectrum phase is 2.5. A E(l) host model has been subtracted. 
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Fig. A.35. SNLS- 04Dlpg684: a SN la supernova at z=0.515. The spectrum phase is -1.3. A Sa-Sb host model has been subtracted. 
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Fig. A.36. SNLS- 04Dlpp683: a SN la supernova at z=0.735. The spectrum phase is 2.2. Best-fit obtained for a model with no host 
galaxy component. 
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Fig. A.37. SNLS- 04Dlqd688: a SN la* supernova at z-Q.161 . The spectrum phase is -0.2. A E(l) host modelhas been subtracted. 
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Fig. A.38. SNLS- 04Dlrh710: a SN la supernova at z=0.436. The spectrum phase is 0.0. Best-fit obtained for a model with no host 
galaxy. 
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Fig. A.39. SNLS- 04Dlrx713: a SN la* supernova at z=0.985. The spectrum phase is 0.9. Best-fit obtained for a model with no 
host galaxy component. 
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Fig. A.40. SNLS- 04Dlsa716: a SN la supernova at z=0.585. The spectrum phase is -2.6. A E(l) host model has been subtracted. 
Note the bad fit in the UV region due to subtracting too blue a host model. 
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Fig. A.41. SNLS- 04Dlsi716: a SN la supernova at z=0.702. The spectrum phase is -1.7. A Sa(9) host model has been subtracted. 
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Fig. A.42. SNLS- 04D2ac395: a SN la supernova at z=0.348. The spectrum phase is 1.3. A Sa-Sb host model has been subtracted. 
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Fig. A.43. SNLS- 04D2al390: a SN la supernova at z=0.836. The spectrum phase is -2.5. A E(2) host model has been subtracted. 




Fig. A.44. SNLS- 04D2an386: a SN la supernova at z=0.62. The spectrum phase is -3.4. Best-fit obtained for a model with no host 
galaxy component. Note that the Ca n feature at rest-frame 3700 A is deeper than in the model. 
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Fig. A.45. SNLS- 04D2an392: a SN la supernova at z=0.62. The spectrum phase is 0.3. A Sc(9) host model has been subtracted. 




Fig. A.46. SNLS- 04D2bt444: a SN la supernova at z=0.220. The spectrum phase is 6.6. A E(6) host model has been subtracted. 




Fig. A.47. SNLS- 04D2ca445: a SN la* supernova at z=0.835. The spectrum phase is 12.0. Best-fit obtained for a model with no 
host galaxy component. 
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Fig. A.48. SNLS- 04D2ca447: a SN la* supernova at z=0.835. The spectrum phase is 13.1. Best-fit obtained for a model with no 
host galaxy component. 
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Fig. A.49. SNLS- 04D2cc446: a SN la* supernova at z=0.838. The spectrum phase is 6.1. A S0(3) host model has been subtracted. 




Fig. A.50. SNLS- 04D2cc448: a SN la* supernova at z=0.838. The spectrum phase is 7.2. A S0(2) host model has been subtracted. 
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Fig. A.51. SNLS- 04D2cf447: a SN la supernova at z=0.369. The spectrum phase is 8.5. Best-fit obtained for a model with no host 
galaxy component. Note that the model is unable to reproduce the feature at rest-frame 4900 A. 
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Fig. A.52. SNLS- 04D2cw445: a SN la* supernova at z=0.568. The spectrum phase is 19.3. A Sb-Sc host model has been sub- 
tracted. The UV part of the model is poorly constrained due to the lack of UV templates at late phases in the version of the SALT2 
training set used in this work. 
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Fig. A.53. SNLS- 04D2fp470: a SN la supernova at z=0.415. The spectrum phase is 1.8. Best-fit obtained for a model with no host 
galaxy component. 
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Fig. A.54. SNLS- 04D2fs470: a SN la supernova at z=0.357. The spectrum phase is 1.7. Best-fit obtained for a model with no host 
galaxy component. 
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Fig. A.55. SNLS- 04D2gc471: a SN la supernova at z=0.521. The spectrum phase is -4.9. A E(l) host model has been subtracted. 
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Fig. A.56. SNLS- 04D2gp475: a SN la supernova at z=0.732. The spectrum phase is 2.7. Best-fit obtained for a model with no host 
galaxy component. 




Fig. A.57. SNLS- 04D2iu497: a SN la* supernova at z=0.70. The spectrum phase is 9.8. Best-fit obtained for a model with no host 
galaxy component. 




Fig. A.58. SNLS- 04D2iu498: a SN la* supernova at z=0.70. The spectrum phase is 10.4. A E(l) host model has been subtracted. 
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Fig. A.59. SNLS- 04D2ja499: a SN la* supernova at z=0.740. The spectrum phase is 9.2. A E(3) host model has been subtracted. 
Note the large feature at rest-frame 4800 A not reproduced by the model. 
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Fig. A.60. SNLS- 04D2ja500; a SN la* supernova at z=0.740. The spectrum phase is 9.8. A E(3) host model has been subtracted. 
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Fig. A.61. SNLS- 04D2mc736: a SN la supernova at z=0.348. The spectrum phase is 6.5. A Sa-Sb host model has been subtracted. 
Three re-calibration parameters have been used to fit this spectrum. 




Fig. A.62. SNLS- 04D4an558: a SN la supernova at z=0.613. The spectrum phase is 8.L Best-fit obtained for a model with no host 
galaxy component. 
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Fig. A.63. SNLS- 04D4bk558: a SN la* supernova at z=0.88. The spectrum phase is 3.1. A E(l) host model has been subtracted. 
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Fig. A.64. SNLS- 04D4bq558: a SN la supernova at z=0.55. The spectrum phase is 2.6. Best-fit obtained for a model with no host 
galaxy component. 
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Fig. A.65. SNLS- 04D4bq562: a SN la supernova at z=0.55. The spectrum phase is 5.1. Best-fit obtained for a model with no host 
galaxy component. 




Fig. A.66. SNLS- 04D4dw566: a SN la* supernova at z= 1.031. The spectrum phase is 2.1. A Sa-Sb host model has been subtracted. 




Fig. A.67. SNLS- 04D4fx589: a SN la supernova at z=0.629. The spectrum phase is -8.4. Best-fit obtained for a model with no 
host galaxy component. 
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Fig. A.68. SNLS- 04D4gz598: a SN la* supernova at z=0.375. The spectrum phase is -5.8. A Sa-Sb host model has been subtracted. 
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Fig. A.69. SNLS- 04D4gz620: a SN la* supernova at z=0.375. The spectrum phase is 10.2. Best-fit obtained for a model with no 
host galaxy component. 
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Fig. A.70. SNLS- 04D4hf598: a SN la* supernova at z=0.94. The spectrum phase is -0.5. Best-fit obtained for a model with no 
host galaxy component. 
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Fig. A.71. SNLS- 04D4ht628: a SN la supernova at z=0.217. The spectrum phase is 6.6. A Sa-Sb host model has been subtracted. 




Fig. A.72. SNLS- 04D4ib629: a SN la supernova at z=0.699. The spectrum phase is 0.7. A SO-Sa host model has been subtracted. 
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Fig. A.73. SNLS- 04D4id628: a SN la supernova at z=0.769. The spectrum phase is 3.0. Best-fit obtained for a model with no host 
galaxy component. 
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Fig. A.74. SNLS- 04D4jr650: a SN la supernova at z=0.47. The spectrum phase is -5.9. Best-fit obtained for a model with no host 
galaxy component. 
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Fig. A.75. SNLS- 04D4ju650: a SN la supernova at z=0.472. The spectrum phase is -2.2. A S0(1) host model has been subtracted. 
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Fig. A.76. SNLS- 04D4jw651: a SN la* supernova at z=0.96L The spectrum phase is 2.2. A SO-Sa host model has been subtracted. 
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Fig. A.77. SNLS- 05Dlcb984: a SN la supernova at z=0.632. The spectrum phase is 4.3. A Sb-Sc host model has been subtracted. 
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Fig. A.78. SNLS- 05Dlck983: a SN la supernova at z=0.617. The spectrum phase is -2.7. A Sb-Sc host model has been subtracted. 




Fig. A.79. SNLS- 05Dlcl999: a SN la* supernova at z=0.83. The spectrum phase is 8.7. A Sd(4) host model has been subtracted. 
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Fig. A.80. SNLS- 05Dldnl002: a SN la supernova at z=0.566. The spectrum phase is -4.7. Best-fit obtained for a model with no 
host galaxy component. 
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Fig.A.81. SNLS- 05Dlhkl072: a SN Ia_pec supernova at z=0.263. The spectrum phase is -5.0. A Sc(5) host model has been 
subtracted. 




Fig. A.82. SNLS- 05Dlhnl073: a SN la supernova at z=0.149. The spectrum phase is -LO. Best-fit obtained for a model with no 
host galaxy component. This is a very reddened SN la due to host interstellar absorption. 




Fig. A.83. SNLS- 05Dlizl069: a SN la* supernova at z=0.86. The spectrum phase is 7.8. A Sd(l) host model has been subtracted. 
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Fig. A.84. SNLS- 05Dlkel094: a SN la supernova at z=0.69. The spectrum phase is 2.2. Best-fit obtained for a model with no host 
galaxy component. 




Fig. A.85. SNLS- 05D2ac738: a SN la supernova at z=0.479. The spectrum phase is 2.0. Best-fit obtained for a model with no host 
galaxy component. Note that the model overestimates the flux in the UV region of the spectrum. 




Fig. A.86. SNLS- 05D2ay746: a SN la* supernova at z=0.92. The spectrum phase is 5.8. A S0(1) host model has been subtracted. 
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Fig. A.87. SNLS- 05D2bt764: a SN la supernova at z=0.68. The spectrum phase is 1.4. Best-fit obtained for a model with no host 
galaxy component. 
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Fig. A.88. SNLS- 05D2bv764: a SN la supernova at z=0.474. The spectrum phase is -0. 1 . A Sb-Sc host model has been subtracted. 
Note that the Ca ii feature at rest-frame 3700 A is narrower than in the model. 




Fig. A.89. SNLS- 05D2bw764: a SN la supernova at z=0.92. The spectrum phase is 1.9. A Sb-Sc host model has been subtracted. 
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Fig. A.90. SNLS- 05D2by775: a SN la* supernova at z=0.891. The spectrum phase is 0.9. A Sc(ll) host model has been subtracted. 
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Fig.A.91. SNLS- 05D2cb775: a SN la* supernova at z=0.427. The spectram phase is -5.8. A Sc(ll) host model has been sub- 
tracted. 
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Fig. A.92. SNLS- 05D2ci795: a SN la supernova at z=0.630. The spectrum phase is 3.6. Best-fit obtained for a model with no host 
galaxy component. 
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Fig. A.93. SNLS- 05D2ci796: a SN la supernova at z=0.630. The spectrum phase is 4.9. Best-fit obtained for a model with no host 
galaxy component. 
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Fig. A.94. SNLS- 05D2ci799: a SN la supernova at z=0.630. The spectrum phase is 6.7. Best-fit obtained for a model with no host 
galaxy component. 
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Fig. A.95. SNLS- 05D2ct799: a SN la* supernova at z=0.734. The spectrum phase is 7.8. A Sd(l 1) host model has been subtracted. 
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Fig. A.96. SNLS- 05D2dt807: a SN la supernova at z=0.574. The spectrum phase is -L7. A E(3) host model has been subtracted. 
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Fig.A.97. SNLS- 05D2dw807: aSN la supernova at z=0.417. The spectrum phase is -5.3. A S0(1) host model has been subtracted. 
Note that the model is unable to reproduce the deep Ca ii feature at rest-frame 3700 A. 




Fig. A.98. SNLS- 05D2dy805: a SN la supernova at z=0.51. The spectrum phase is LL A Sb-Sc host model has been subtracted. 
Note that the Ca ii feature at rest-frame 3700 A is narrower than in the model. 
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Fig.A.99. SNLS- 05D2eb805: a SN la supernova at z=0.534. The spectrum phase is -4.7. A S0(1) host model has been subtracted. 
Note that the model is unable to reproduce the deep Ca ii feature at rest-frame 3700 A. 
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Fig. A.IOO. SNLS- 05D2ec805; a SN la supernova at z=0.640. The spectrum phase is 2.7. Best-fit obtained for a model with no 
host galaxy component. 
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Fig. A.lOl. SNLS- 05D2ei833: a SN la* supernova at z=0.366. The spectrum phase is 16.9. A S0(4) host model has been sub- 
tracted. 
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Fig. A.lOl. SNLS- 05D2fq827: a SN la supernova at z=0.733. The spectrum phase is 2.1. Best-fit obtained for a model with no 
host galaxy component. 
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Fig. A.103. SNLS- 05D2he833: a SN la supernova at z=0.608. The spectrum phase is 3.0. A E(l) host model has been subtracted. 
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Fig. A.104. SNLS- 05D2ie829: a SN la supernova at z=0.348. The spectrum phase is -8.9. A Sc(4) host model has been subtracted. 
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Fig. A.105. SNLS- 05D2nnl071: a SN la supernova at z=0.87. The spectrum phase is -1.3. A Sa(8) host model has been subtracted. 
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Fig. A.106. SNLS- 05D4af914: a SN la supernova at z=0.499. The spectrum phase is 1 1.3. Best-fit obtained for a model with no 
host galaxy component. 
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Fig. A.107. SNLS- 05D4ag922: a SN la* supernova at z=0.64. The spectrum phase is 14.0. A Sb-Sc host model has been subtracted. 
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Fig. A.108. SNLS- 05D4ay914: a SN la* supernova at z=0.408. The spectrum phase is 8.2. A Sa-Sb host model has been subtracted. 
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Fig. A.109. SNLS- 05D4ay9 16: a SN la* supernova at z=0.408. The spectrum phase is 9.6. A Sa-Sb host model has been subtracted. 
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Fig. A.llO. SNLS- 05D4be918: a SN la supernova at z=0.537. The spectrum phase is 3.9. A E(l) host model has been subtracted. 
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Fig. A.m. SNLS- 05D4bi918: aSN la supernova at z=0.775. The spectrum phase is -1.5. A Sc(l) host model has been subtracted. 
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Fig. A.112. SNLS- 05D4bj921: a SN la supernova at z=0.701. The spectrum phase is -2.0. A Sc(9) host model has been subtracted. 
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Fig.A.113. SNLS- 05D4cn952: aSN la supernova at z=0.763. The spectrum phase is 4.7. A S0(1) host model has been subtracted. 




Fig. A.114. SNLS- 05D4cq946: a SN la supernova at z=0.70L The spectrum phase is -O.L Best-fit obtained for a model with no 
host galaxy component. 
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Fig. A.115. SNLS- 05D4cs946: a SN la supernova at z-0.19. The spectrum phase is -1.7. A Sc(12) host model has been subtracted. 
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Fig. A.116. SNLS- 05D4cw950: a SN la supernova at z=0.375. The spectrum phase is 7.0. A S0(6) host model has been subtracted. 
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Fig. A.117. SNLS- 05D4dw983: a SN la supernova at z=0.855. The spectrum phase is 4.7. A Sc(12) host model has been subtracted. 




Fig. A.118. SNLS- 05D4ef997: a SN la* supernova at z=0.605. The spectrum phase is 3.5. A E(3) host model has been subtracted. 



Balland et al.: The ESOA'LT 3rd year Type la supernova data set from the SNLS, Online Material p 62 

Rest-frame wavelength (A) 

„ 2500 3000 3500 4000 4500 5000 5500 
2.5l I 1 1 I I I 1 I I I I 1 1 1 1 I 1 1 I I I 1 I I I I 1 1 1 1 I 1 1 




Observed Wavelength (A) 



Fig. A.119. SNLS- 05D4ejl000: a SN la supernova at z=0.585. The spectrum phase is 7.6. Best-fit obtained for a model with no 
host galaxy component. 




Fig. A.120. SNLS- 05D4ek999: a SN la supernova at z=0.536. The spectrum phase is 2.1. Best-fit obtained for a model with no 
host galaxy component. 
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Fig. A.121. SNLS- 05D4ev985: a SN la supernova at z=0.722. The spectrum phase is -3.6. Best-fit obtained for a model with no 
host galaxy component. 




Fig. A.122. SNLS- 05D4fel000: a SN la* supernova at z=0.984. The spectrum phase is -2.0. A E(3) host model has been subtracted. 
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Fig. A.123. SNLS- 05D4fF1003: a SN la* supernova at z=0.402. The spectrum phase is 5. 1 . A Sc(6) host model has been subtracted. 
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Fig. A.124. SNLS- 05D4fgl002: a SN la supernova at z=0.839. The spectrum phase is -0.3. A S0(5) host model has been subtracted. 
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Fig. A.125. SNLS- 06Dlabl 100: a SN la supernova at z=0.182. The spectrum phase is -4.4. A E(7) host model has been subtracted. 
Note that the model is unable to reproduce the large feature at rest-frame 4500 A. 
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Fig. A.126. SNLS- 06D2bkl 135; a SN la supernova at z=0.499. The spectrum phase is 0.9. A Sb-Sc host model has been subtracted. 



Balland et al.: The ESOA'LT 3rd year Type la supernova data set from the SNLS, Online Material p 66 




Fig. A.127. SNLS- 06D2cal 153: a SN la supernova at z=0.531. The spectrum phase is 0.0. A Sd(5) host model has been subtracted. 




Fig. A.128. SNLS- 06D2cal 154: a SN la supernova at z=0.531. The spectrum phase is 0.7. A Sd(3) host model has been subtracted. 
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Fig. A.129. SNLS- 06D2cbl 155: a SN la* supernova at z=1.00. The spectrum phase is 3.0. Best-fit obtained for a model with no 
host galaxy component. 
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Fig. A.130. SNLS- 06D2ccl 153: a SN la supernova at z=0.532. The spectrum phase is 3.3. Best-fit obtained for a model with no 
host galaxy component. 
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Fig. A.131. SNLS- 06D2cdl 157: a SN la supernova at z=0.930. Spectrum obtained with the 3001 Grism on FORS2. The spectrum 
phase is 4.2. Best-fit obtained for a model with no host galaxy component. 
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Fig. A.132. SNLS- 06D2cel 156: a SN la supernova at z=0.82. Spectrum obtained with the 3001 Grism on FORS2. The spectrum 
phase is 0.0. Best-fit obtained for a model with no host galaxy component. 
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Fig.A.133. SNLS- 06D2ckll62: a SN la* supernova at z=0.552. The spectrum phase is 7.7. A Sa-Sb host model has been 
subtracted. 
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Fig. A.134. SNLS- 06D2gal241: a SN la* supernova at z=0. 84. The spectrum phase is 5.4. A Sb-Sc host model has been subtracted. 
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Fig. A.135. SNLS- 06D4cel295: a SN la supernova at z=0.85. The spectrum phase is 2.7. Best-fit obtained for a model with no 
host galaxy component. 
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Fig. A.136. SNLS- 06D4cel296: a SN la supernova at z=0.85. The spectrum phase is 3.3. Best-fit obtained for a model with no 
host galaxy component. 
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Fig. A.137. SNLS- 06D4cll297: a SN la supernova at z=1.00. The spectrum phase is -2.5. A E(l) host model has been subtracted. 
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Fig. A.138. SNLS- 06D4col306: a SN la supernova at z=0.437. The spectrum phase is 3.5. A E-SO host model has been subtracted. 
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Fig. A.139. SNLS- 06D4cql306: a SN la supernova at z=0.41 1. The spectrum phase is -1.4. Best-fit obtained for a model with no 
host galaxy component. 



